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Abstract 

We study nonlinear systems of hyperbolic (in a wider sense) PDE's in R'* describing 
wave propagation in dispersive nonlinear media such as, for example, electromagnetic 
waves in nonlinear photonic crystals. The initial data is assumed to be a finite sum 
of wavepackets referred to as a multi-wavepacket. The wavepackets and the medium 
nonlinearity are characterized by two principal small parameters l3 and g where: (i) i 
is a factor describing spatial extension of involved wavepackets; (ii) is a factor de- 
scribing the relative magnitude of the linear part of the evolution equation compared to 
its nonlinearity. A key element in our approach is a proper definition of a wavepacket. 
Remarkably, the introduced definition has a flexibility sufficient for a wavepacket to 
preserve its deflning properties under a general nonlinear evolution for long times. In 
particular, the corresponding wave vectors and the band numbers of involved wavepack- 
ets are "conserved quantities". We also prove that the evolution of a multi-wavepacket 
is described with high accuracy by a properly constructed system of envelope equations 
with a universal nonlinearity. The universal nonlinearity is obtained by a time averag- 
ing applied to the original nonlinearity, in simpler cases the averaged system turns into 
a system of Nonlinear Schrodinger equations. 

1 Introduction 

The underlying physical subject of this work is propagation of a multi-wavepacket (a finite 
system of wavepackets) in a spatially dispersive and nonlinear medium, and we are particu- 
larly interested in electromagnetic waves propagation in nonlinear photonic crystals, see [55] . 
|56j . [58] . [4j-[7] and references therein, with the nonlinear optics constitutive relations, [12], 
[T51 Sections 1,2], [32], [33 • The mathematical subject of interest is the following general 
nonlinear evolutionary system 

drV = --L (-iV) U + F (U) , U (r, t)1_^ = h (r) , r e (1) 
g 

where (i) U = U (r, r), r G M"^, U G C^"^ is a 2 J dimensional vector; (ii) L (— iV) is a lin- 
ear self-adjoint differential (pseudodifferential) operator with constant coefficients with the 
symbol L (k), which is a Hermitian 2 J x 2 J matrix; (iii) F is a polynomial nonlinearity such 
that F (0) = 0, F' (0) = and F (U) is translation- invariant, i.e. if TJJ (r) = U (r a) for 
aeR'^ then F (TJJ) = T^F (U); (iv) h = h (r) is assumed to be the sum of a finite number 
of wavepackets h/, I — 1, . . . , iV; (v) > is a small parameter. In the case of nonlinear 
photonic crystals the components of the vector field U (r) are the modal amplitudes of the 
electromagnetic field and the nonlinearity F (U) is constructed from the nonlinear medium 
polarization in the adiabatic approximation, [15, Sections 2.4.2]. The systems of the form 
([T]) also describe as a particular case well-known equations, namely: complexification of the 
Nonlinear Schrodinger equation; coupled envelope equations which arise in nonlinear bire- 
fringent optical media, [H] Section 2i]; nonlinear Klein-Gordon and Sine-Gordon equations 
[611, Section 14.1], [331 Section 5.8.3], [331 Section 9.6]. Such equations appear in a number 



of physical problems: elementary particles; dislocations in crystals; propagation of Bloch's 
domain walls in the theory of ferromagnetism; self-induced transparency in nonlinear optics; 
the propagation of magnetic flux quanta in long Josephson transmission lines. Significance 
and importance of wavepacket solutions from the both physical and mathematical points of 
view is discussed in [3-0, lU Section 2], [S^, fSSj. 

There are numerous problems involving small parameters only in the initial data which 
can be reduced to the form ([!]), for instance, problems with high frequency initial data or 
small initial data with consequent evolution on long time intervals (see Section 3 for details) . 

We study the nonlinear evolution equation ([1]) on a finite time interval 

< T < r», where r» > is a fixed number. (2) 

The time r, may depend on the L°° norm of the initial data h but, importantly, does 
not depend on g. We consider classes of initial data such that wave evolution governed by 
([IP is significantly nonlinear on time interval [0, t^] and the effect of the nonlinearity F (U) 
does not vanish as g — > 0. 

Since the both linear operator L (— iV) and the nonlinearity F (U) are translation in- 
variant, it is natural and convenient to recast the evolution equation ([1]) by applying to it 
the Fourier transform with respect to the space variables r, namely 

9,U(k) = --L(k)U(k)-fFfu)(k), U(k) =h(k), (3) 

g \ J T=o 

where U (k) is the Fourier transform of U (r), i.e. 

U(k)=/" U(r)e-''" ''dr, U(r) = (27r)"'' /" U (k) e''" ''dr, where r, k e M'', (4) 

and F is the Fourier form of the nonlinear operator F (U) involving convolutions. 

The nonlinear evolution equations ([T]), Q are commonly interpreted as describing wave 
propagation in a nonlinear medium. We assume that the linear part L (k) is a 2 J x 2 J 
Hermitian matrix with eigenvalues iOn,Q (k) and eigenvectors g„.^ (k) satisfying 

L (k) ^r^x (k) = w«,c (k) g„,c (k) , C = ±, ^n.+ (k) > 0, c^„,_ (k) < 0, n - 1, . . . , J, (5) 

where uJn,Q (k) are real-valued, continuous for all non-singular k functions, and vectors 
gn,C (k) G C^'^ have unit length in the standard Euclidean norm. The functions LOn.c, (k), 
n — 1, . . . , J, are called dispersion relations between the frequency oj and the wavevector k 
with n being the hand number. We assume that the eigenvalues are naturally ordered by 

WX+ (k) > . . . > (k) > > LJi.- (k) > . . . > uj.- (k) , (6) 

and for almost every k (with respect to the standard Lebesgue measure) the eigenvalues 
are distinct and, consequently, the above inequalities become strict. Importantly, we also 
assume the following diagonal symmetry condition 

w„ -c (^k) = -uJnX (k) , C = ±, = 1, • • ■ , (7) 

which is naturally present in many physical problems (see also Remark [TJ] below) , and is a 
fundamental condition imposed on the matrix L (k). In addition to that in many examples 
we also have 

g„^^ (k) ~ Zn-c (^k) , where z* is complex conjugate to z. (8) 
Very often we will use the following abbreviation 

W„,+ (k) = UJn (k) . (9) 
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From we obtain 



(k) = -c^„ (-k) , uj„^c (k) = C^n (Ck) , C = ±- (10) 

We also will often use the orthogonal projection n„_^ (k) in C^'' onto the complex line defined 
by the eigenvector g„_^ (k), namely 

Unx (k) u (k) = t2„,c (k) gn,c (k) = UnX (k) , n = 1, . . . , J, C = ±. (11) 

As it is indicated by the title of this paper we study the nonlinear problem ([T]) for initial 
data h in the form of a properly defined wavepacket or, more generally, a sum of wavepackets 
which we refer to as multi-wavepacket. The simplest example of a wavepacket w is provided 
by the following formula 

w (r, (3) = {(3r) e"'---g„,+ (k,) , r G (12) 

where k* S R'^ is a wavepacket wave vector, n is hand number, and /3 > is a small pa- 
rameter. We refer to the pair (n,k,) in (|12p as wavepacket nk-pair. Observe that the 
space extension of the wavepacket w (r,/3) is proportional to /3~^ and it is large for small 
(3. Notice also that if /3 — > the wavepacket w (r, (3) as in tends, up to a constant 
factor, to the elementary eigenmode e''^* '"g„^i^ (k*) of the operator L (— iV) with the corre- 
sponding eigenvalue Wn,c (k*). We refer to wavepackets of the simple form p2p as simple 
wavepackets to underline the very special way the parameter (3 enters its representation. 
The function <i>^ (r), which we call wavepacket envelope, describes its shape and it can be 
any scalar complex-valued regular enough function, for example a function from Schwartz 
space. Importantly, as /3 — s- the L°° norm of a wavepacket US\) remains constant, and, 
hence, nonlinear effects in (QP remain strong. 

Evolution of wavepackets in problems which can be reduced to the form ^ were studied 
for a variety of equations in numerous physical and mathematical papers, mostly by asymp- 
totic expansions with respect to a single small parameter similar to /3, see [TU], [T3], [TB], 
ppj . P] . [3D], [35], HZ], [SD], [5T] and references therein. We are interested in general 
properties of evolutionary systems of the form ([T]) with wavepacket initial data which hold 
for a wide class of nonlinearities and all values of the space dimensions d of the number 2 J of 
the system components. Our approach is not based on asymptotic expansions but involves 
the two small parameters (3 and g with mild constraints on their relative smallness. The 
constraints can be expressed either in the form of certain inequalities or equalities, and a 
possible simple form of such a constraint can be a power law 

f3 — Cq'^ where C > and >t > Q are arbitrary constants. (13) 

Of course, general features of wavepacket evolution are independent of particular values 
of the constant C In addition to that, some fundamental properties such as wavepacket 
invariance, are also totally independent on particular choice of the values of x in p3p . 
whereas other properties are independent of x as it varies in certain intervals. For for 
instance, dispersion effects are dominant for x < 1/2, whereas the wavepacket superposition 
principle of T holds for x < 1. 

The qualitative picture of wavepacket evolution dependence on small (3 and g is as follows. 
The parameter /3 enters the problem ^ through the multi-wavepacket initial data h (r, (3) 
whereas g enters it through the factor i before the linear part. Evidently the factor ^ 
determines the relative magnitude of the linear part compared to the nonlinearity and since 
i is large, one expects the linear part to provide an important input into solutions properties. 
This input includes, in particular, key role of eigenmodes and eigenfreguencies (dispersion 
relations) in expressing the nonlinear evolution. Importantly, in many cases of interest 
though i is large, nonlinear phenomena are significant and this is the case when (3 < C g^^^ . 

More precisely, if /3 < Cg^l"^ then, as in the case of finite-dimensional nonlinear ODE 
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evolutionary systems, the large values of ^ lead to a well defined solution factorization into 
the fast (high frequency) and the slow (low frequency) components. The interplay between 
the fast and slow components is also similar to the ODE case, namely, the nonlinear evolution 
is associated primarily with the slow component governed by a nonlinear equation obtained 
from the original one by a certain canonical time averaging procedure. Our further analysis of 
the above mentioned interplay shows the following. Firstly, the linear superposition principle 
holds, [7] , that is if < 1 is as in (|13p and the initial data is a sum of generic wavepackets then 
the solution is the sum of the solutions for single involved wavepackets with precision -^Stt 
with arbitrary small e. Secondly, properly defined wavepackets and their linear combinations 
are preserved under the nonlinear evolution ([T]), which is a subject of this paper. 

In the light of the above discussion we introduce the slow variable u (k, t) by the formula 

U(k,T) =e-TL('')u(k,r), (14) 
and recast the equation ([3]) for it as follows 

dru = ef (e^^u) , u|^^o = h. (15) 

Then we obtain an integral form of (fT5|l by integrating it with respect to r: 

u = J- (u) + h, T{u)=T {g) (u) = e'-^^F (e^^u (r')) dr' (16) 

with explicitly defined nonlinear polynomial integral operator J-{q), which depends on the 
parameter g. This operator is bounded uniformly with respect to g in the Banach space 
-B = C ([0, r,] , L^) of functions v (k, r), < r < r*, with the norm 

||v(k,T)||^ = ||v(k,r)I|p([o..]Li)= sup / |v(k,T)|dk, (17) 

0<r<r. jRd 

where is the Lebesgue space of functions v (k) with the standard norm 

||v(.)ILi- / |v(k)|dk. (18) 



Sometimes we use more general weighted spaces L^'° with the norm 

l|v|L.. = / (l + |k|)"|v(k)|dk, a>0. (19) 

A rather elementary existence and uniqueness theorem fTheoremE^ implies that for a small 
and, importantly, independent of g constant > this equation has a unique solution 

Hr)^Q [q) , h) (t) , t e [0, r,] , u e ([0, r,] , L^) , (20) 



where Q denotes the solution operator for the equation p6p . the operator depends on oper- 
ator T (g), which itself depends on the parameter g. If u (k, r) is a solution to the equation 
((T6)) we call the function U (r, r) defined by ^ an F-solution to the equation ([T]). We 
denote by the space of functions V (r) such that their Fourier transform V (k) belongs 



to L^, and define ||V||j^i 



. Since 



V 

||V||^„ < (271)^" V and C i°° (21) 
i^-sohitions to 11]) belong to ([0,t*] ,2^) C ([0, t*] , L°°). 
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We would like to define wavepackets in a form which explicitly allows them to be real 
valued. This is accomplished based on the symmetry ([7]) of the dispersion relations by 
introduction of a doublet wavepacket 

w (r, P) = iPr) e*---g„,+ (k,) + (/3r) e-"'---g„,_ (-k,) . (22) 

Such a wavepacket is real if $_ (r) ,g„,- k») is complex conjugate to $+ (r), g„^+ (k,), 
i.e. if 

$ _ (r) = (r) , g„,+ (k,) = g„,_ (-k,)* . (23) 

Considering wavepackets with nfc-pair (n,k,) we usually mean doublet ones as in (|22p . but 
sometimes or $_ may be zero producing (jl2[) . 

To identify characteristic properties of a wavepacket suitable for our needs, let us look 
at the Fourier transform w (k, /3) of an elementary wavepacket w (r, /3) defined by (HH), that 
is 

w (k, P) = /3-'*<i. (k - K)) g„,c (k,) . (24) 

We call such w (k, f3) wavepacket too, obviously it possesses the following properties: (i) its 
norm is bounded (in fact, constant), uniformly in /3 — > 0; (ii) for every e > the value 
w (k, /3) — > for every k outside a /3^~'^-neighborhood of k», and the convergence is faster 
than any power of /3 if <i> is a Schwartz function. To explicitly interpret the last property we 
introduce a cutoff function [rf) 

* (J7) = 1 for \ri\ < 1, * (J7) = for \ri\ > 1, (25) 

together with its shifted/rescaled modification 

^ik-K)^-^{k;K,l3^-') =*(/3-(i^^)(k-kO) . (26) 

If in an elementary wavepacket w (r, /?) defined by (I24p (f>^ (r) is a Schwartz function then 

||(l-*(-,k„/3i-^))w(.,/3)|| <C,,,r, 0</3< 1, 

which holds for arbitrarily small e > and arbitrarily large s > 0. Based on the above 
discussion we give the following definition of a wavepacket which is a minor variation of [71 
Definiton 8]. 

Definition 1 (single-band wavepacket) Let < e < 1 be a fixed number. For a given 
band number n G {1, . . . , J} and a wavevector k* G R'^ a function h(/3,k) is called a 
wavepacket with nfc-pair (n,k*) and the degree of regularity s > if there exists such 
/3g > that for (3 < Pq the following conditions are satisfied: (i) h (/3, k) is L^-hounded 
uniformly in fi, i.e. 

h (/3, •) ^^<C, < /3 < /?o for some C > 0; (27) 

(ii) h (/?, k) has the following structure: 

h {(3, k) = h_ (/3, k) + h+ (/?, k) + < /3 < /?o, where (28) 

he (/?, k) = * (k, Ck, , (3^-') n„,c (k) he (/3, k) , C = ±, (29) 
with (•, ^k*, /3"'^~'^) defined by i26\) and Dh satisfying the following tail estimate.' 



< C'(3% 0<(3 <(3„ for some C > 0. (30) 



The inverse Fourier transform h (/3, r) of a wavepacket h (/9, k) is also called a wavepacket. 
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Point (ii) of the above Definition means that the wavepacket h (/3, k) is composed of 
two functions h^; (/3,k), ^ = ±, which take values in n-th band eigenspace of L (k) and are 
localized near Ck*, where (n, k*) is the nfc-pair of the wavepacket. The number /3q usually 
is small and may depend on a wavepacket. 

Evidently, if a wavepacket has the degree of regularity s, it also has a smaller degree of 
regularity s' < s with the same e. Observe that the degree of regularity s is related to the 
smoothness of (r) in p2|) so that the higher is the smoothness the higher I can be taken. 
Namely, if £ L^'" then one can take any I < a, see Lemma [5^ below. For example, if 
in the elementary wavepacket w (r, f3) defined by $^ (r) is a Schwartz function then it 
has arbitrarily large degree of regularity. 

It turns out remarkably that wavepackets satisfying Definition [1] preserve their defining 
properties under nonlinear evolution. It is remarkable, in particular, since it is well-known 
that determination of classes of solutions which preserve their form under generic nonlinear 
evolution usually leads to infinite expansions, such as multi-scale expansions, power expan- 
sions, modal expansions etc with serious difficulties in establishing the convergence. Such 
expansions often are formally invariant, but they involve infinitely many rather complex 
terms and establishing the convergence is a very hard problem indeed if there is any con- 
vergence at all. Our Definition [T] of a wavepacket involves only a finite number of terms 
and its invariance is provided by the flexible tail term Df^. We also find remarkable the 
very simplicity of the definition which nevetherless allows for a sufficiently detailed analysis 
of the dynamics, including, in particular, rigorously justified NLS-type approximations of 
wavepacket dynamics presented in the following sections. 

Our special interest is in waves that are finite sums of wavepackets and we refer to them 
as multi-wavepackets. 

Definition 2 (multi-wavepacket) Let S be a set of nk-pairs: 

S^{{ni,Ki), ; = l,...,Ar}cS = {l,...,J}xR^ {ni,Ki)^{ni,,Ki,) for I ^ I', (31) 

and N — \S\ he their number. Let Ks be a set consisting of all different wavevectors k,; 
involved in S with \Ks\ < N being the number of its elements. Ks is called wavepacket 
fc-spectrum and without loss of genericity we assume the indexing of elements in S to be 
.such that 

= {k„,z = l,...,|Xs|}, *.e. k^iforl<i< \Ks\. (32) 

A function h (/3) = h (/3, k) is called a multi-wavepacket with nfc-spectrum S if it is a finite 
sum of wavepackets, namely 

N 

h (/?, k) = ^ k) , < /? < ^0 for some (3q > 0, (33) 
1=1 

where h.i, I = I, . . . , N , is a wavepacket with nk-pair (k*;, ni) & S as in Definition\^ 

Note that if h (/?, k) is a wavepacket then h (/?, k) + O {(3^) is a wavepacket as well with 
the same nfc-spectrum, and the same is true for multi-wavepackets. Hence, we can introduce 
multi-wavepackets equivalence relation " ~" of the degree s by 

hi (/3, k) ~ ha (/3, k) if hi k) - ha (/3, k) < for some constant C > 0. (34) 

Let us turn now to the abstract nonlinear problem (|16p where (i) J- = J- {g) depends 
on Q and (ii) the initial data h = h (/3) is a multi-wavepacket depending on /3. We would 
like to state our first theorem on multi-wavepacket preservation under the evolution (|16|) for 
/3, g — * 0, which holds, as it turns out, provided its n/c -spectrum S satisfies certain natural 
condition called resonance invariance. This condition is intimately related to the so-called 
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phase and frequency matching conditions for stronger nonlinear interactions, and its concise 
formulation is as follows. We define for given dispersion relations {a;„ (k)} and any finite 
set S" C {1, . . . , J} X R'' another finite set 7^ (5') C {1, . . . , J} x R'^ where 7^ is a certain 
algebraic operation described in Definition [TS] below. It turns out that for any S always 
S C TZ (S) but if, in fact, TZ (S) = S we call S resonance invariant. The condition of 
resonance invariance is instrumental for the multi-wavepacket preservation, and there are 
examples showing that if it fails, i.e. TZ (S) ^ S, the wavepacket preservation does not 
hold. Importantly, the resonance invariance TZ (S) — S allows resonances inside the multi- 
wavepacket, that includes, in particular, resonances associated with the second and the third 
harmonic generations, resonant four-wave interaction etc. 

Theorem 3 (multi-wavepacket preservation) Suppose that the nonlinear evolution is 
governed by il6]} and the initial data h = h(/?,k) is a multi-wavepacket with nk-spectrum 
S and the regularity degree s, and assume S to he resonance invariant (see Definition \18\ 
below). Let dependence between parametrs g and [3 be given be any function g = p {(3) 
satisfying 

< p (/3) < CP'', for some constant C > 0, (35) 

and let us set g — p{l3). Then the solution u(t,/3) — Q ^JF (p (/3)) , h (/?)^ (t) to \10j) for 
any t S [0,r*] is a multi-wavepacket with nk-spectrum S and the regularity degree s, i.e. 

u (r, 13; k) = i 1^'^ ^) ' "^^^^^ wavepacket with nk-pair {ni, k*;) £ S. (36) 

The time interval length > depends only on L^-norms of hi (/3,k) and N. The presen- 
tation i36]) is unique up to the equivalence (3^. 



The above statement can be interpreted as follows. Modes in n/c-spectrum S are always 
resonance coupled with modes in TZ (S) through the nonlinear interactions, but if TZ {S) — S 
then (i) all resonance interactions occur inside 5* and (ii) only small vicinity of S is involved 
in nonlinear interactions leading to the multi-wavepacket preservation. 

Many nonlinear evolution problems with small initial data can be readily reduced by 
elementary rescaling to the system ([T]) with a large parameter ^ before its linear part. For 
example, suppose that F (V) is a homogeneous nonlinearity of degree m (m = 3 for cubic 
one) and that the nonlinear evolution is governed by 

atV = -iL(-iV)V + F(V), V(r,t)|,^o = £''^^™"'^h(r), reM^ (37) 

considered for small g on the large time interval < t < ^ with a fixed > 0. Then the 
following simple change of variables 

V(t) = ei/("-i)U(T), T = tg (38) 

transforms the problem (I37p into equivalent problem In this case the inequality psp 
describes a constraint between spatial extension i and the amplitude factor = 

p (Z?)^^*-'""^' of the initial data. Observe that the equation ([57]) does not have any small 
parameters and the both small parameters g and /3 enter the problem through its initial 
data. Theorem [3] can be restated for the problem (f37|) as follows. 

Corollary 4 (multi-wavepacket preservation) Let V (r, t) be a solution to the nonlin- 
ear system p(/3) is as in I135\) and we set g ~ pif3)- Then if the initial data is such 
that gi~i/(™-i)V (k, 0) = h (k) is a multi-wavepacket then g^-'^/^'^^-'^^V (k, t) remains to be 
a multi-wavepacket with the same nk-spectrum and the degree of regularity for all times 
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The statements of Theorems [3] and Corollary Indirectly follow from the following general 
theorem which makes no assumptions on the relations between f3, g 0. 

Theorem 5 (multi-wavepacket approximation) Let the initial data h in the integral 
equation I116\) be a multi-wavepacket h k) with nk-spectrum S as in iSl]). the regularity 
degree s and with the parameter e > as in Definition [7J Assume that S is resonance 
invariant in the sense of Definition \18\ below. Let the cutoff function ^'(k, k,) and the 
eigenvector projectors Tln.± (k) be defined by 1126]) and Ul]) respectively. For a solution u of 
I116\) we set 

Ul (r, /3; k) = * (k, Ck,;) n„,,c (k)] u (r, k) , ^ = 1, . . . , TV. (39) 

Then every such U; (k;r,/3) is a wavepacket and 



sup 

0<r<T* 



u(T,/3;k)-^^^^u, (T,/3;k) 



<Cig + C2P'' (40) 



where the constant Ci does not depend on e, s and P, and the constant C2 does not depend 
on p. 

It is interesting to note that the statement of the Theorem [5] can be extended to the 
special limit case P — 0, k*; = 0. In this case the initial data of ^ are constants in r 
and we can consider solutions U ([T]) which do not depend on r. Then VU — 0, the linear 
operator L (— iV) reduces to the multiplication by a matrix Lq = L (0) and the system ([1} 
turns into a system of ordinary differential equations (ODE). Notice that (i) the structure of 
the eigenvalues ([7]) implies that the linear part is time-reversible; (ii) the nonlinear part can 
be an arbitrary polynomial. The extension of TheoremOto this case (see Theorem [TT]) reads 
that in a generic, non-resonant situation if initial data are bounded and a set of eigenmodes 
of the matrix Lq is excited at r = then in the course of evolution on a time interval [0, r*] 
where depends on magnitude of initial data (i) all remaining modes remain unexcited 
with accuracy proportional to g, and (ii) only the originally excited modes can significantly 
evolve with this level of accuracy. For finite-dimensional systems governed by ODE's such a 
statement can be derived from the classical time-averaging principle and the time-averaged 
equations remain nonlinear. For infinitely-dimensional systems governed by PDE and with 
the linear operator having a continuous spectrum, as in Theorem [Sj the analysis is more 
complex but the time-averaging still plays important role yielding an accurate approximation 
governed by a certain universal nonlinear PDE. 

We would like to point out also that though Theorem [3] is a simple corollary of more 
general Theorem [5l it is important that the statement 114 Oj l can be formulated as multi- 
wavepacket invariance. That, in particular, allows to take values u(t<,) as new wavepacket 
initial data for ([1]) and extend the wavepacket invariance of a solution to the next time 
interval r, < r < t,i. This observation allows to extend the wavepacket invariance to 
larger values of r (up to blow-up time or infinity) if some additional information about 
solutions with wavepacket initial data is available. In particular, the following theorem 
holds. 



Theorem 6 Assume that all conditions of Theorem\^are satisfied and, in addition to that, 
solutions u (r) of with the multi-wavepacket initial data h (/?) exist on a maximal in- 
terval < T < To such that ||u||pQQ ^i-j < R (ti) for any ti < tq < 00 where R (ti) does 

not depend on P,g. Then the solution u (r, /?) = Q (^J- (p {(3)) ,h.{f3)j (r) to il6\) for any 

T < tq is a multi-wavepacket with nk-spectrum S and the regularity degree s, that is f36\) 
holds. 
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Note that the wavepacket form of solutions can be used to obtain long-time estimates 
of solutions. Namely, very often behavior of every single wavepacket is well approximated 
by its own nonlinear Schrodinger equation (NLS), see [TT], [M], [H], [13], [30], [31], [IZ], 
[50] . [5T] . [53] and references therein, see also Section 6. Many features of the dynamics 
governed by NLS-type equations are well- understood, see [TJ, [Hj, [32, [IH], [SI], [52] and 
references therein. These results can be used to obtain long-time estimates for every single 
wavepacket (as, for example, in [31 ) and, with the help of the superposition principle, for 
the multiwavepacket solution. 

The wavepacket representation (I3S)) from Theorem |3] can be used for more detailed 
analysis of dynamics of wavepackets u; (r, (3) and interaction between them. The following 
theorem illustrates that by describing wavepacket interaction based on a system with a 
weakly universal nonlinearity similar to so-called coupled modes systems or NLS. 

Theorem 7 (NLS-type approximation) Let the conditions of Theorem hold and, in 
addition to that, the initial data hi (k) are of the form h; = h(,_|- + h;._ + Di where 

h,,c(k)=/3-'i?Lc(r'(k-Ck*i))gn„c(k) for \k-Ki\<(3'-\ C = ±, 

Di satisfies i30\) . and every function i?;.,^ {rj), which may depend on (3, is defined for all r] 
and is bounded in L^'"' with a > f uniformly in (3. Then one can write a nonlinear system of 
differential equations for 2N scalar envelope functions zi^q (t, r) with the initial data Hi^q, 
a linear part of the system has order fi < 3 and the nonlinearity is weakly universal as in 
Ii238]) and has order v <1. Let zi^q (r, k), Z = 1, N , be the Fourier transform of a solution 
to this system Then there exist /3q > and a constant C which does not depend on [3, g such 
that for (3 < (3q the solution u of with initial data h can be approximated as follows: 



N 

" <C 



1=1 



E 



Q 



(41) 



The above-mentioned system with a weakly universal nonlinearity is constructed based 
on the equation ([T]) and nfc-spectrum 5* with the help of time averaging l|70)) described below. 
Note that in the simplest case when ij. = 2,v = 0, N= \ (and J is arbitrary) the resulting 
system with a universal nonlinearity is equivalent to classical Nonlinear Schrodinger equa- 
tion (NLS). If iV = 2 and k*i = — k*2 we obtain well-known coupled modes system for 
counterpropagating waves. This theorem applied to particular systems implies approxima- 
tion theorems similar to results of (i) [30] , [53] , [6] , [23] on NLS approximation; (ii) [6] , [24] , 
[47] . [52] on coupled mode approximation; (iii) [54] on three- wave approximations. Note also 
that (|4T]) implies that if p = with 1 < x:' < 2, then the both the first order hyperbolic 
equations (/i = 1, = 0) and the second-order NLS {pL = 2, v = Q) provide an approximation 

for a solution u of (flB]) . but NLS provides a better approximation O {j3^'^"~'^^^ compared with 

O {j3^^^^'^^^^ ^ for first order hyperbolic equations. 

Observe that in the form (|22p for a simple wavepacket we require gn,± (k*) to be an 
eigenvector of the Hermitian matrix L (k*), and one can wonder if gn,± (k*) can be replaced 
with an arbitrary pair of vectors g± in the case J > 1. The answer is affirmative, since one 
can always expand any g with respect to the basis gn.± (k) using n„.± (k), but the result 
will be a multi-wavepacket with up to 2J components rather than a single wavepacket. 

The rest of the paper is organized as follows. In the next section we illustrate important 
points of parameter dependence and wavepacket preservation based on examples. In Section 
3 wc formulate conditions of wavepacket preservation including the key resonance invariance 
condition. In sections 4 we provide examples of different forms of equations and systems 
which involve small or large parameters and can be written in the form of (|T]) after a rescaling. 
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In section 5 we introduce and discuss integrated modal forms of the evolution equation. In 
Section 6 we introduce and study the wavepacket interaction system in its relation to the 
original system. In Section 7 we approximate the wavepacket interaction system by a certain 
minimal wavepacket interaction system which in simplest cases turns into the NLS or the 
coupled modes system. 

2 Preliminary discussion and examples 

Observe that the multi-wavepacket preservation as described in Theorems [3][7] states in 
different forms that (i) its modal composition is essentially preserved; (ii) its nfc-spectrum 
(the set of nA;-pairs {k*;,ni}) remains the same at all times; (iii) no new modes are excited 
with a good accuracy as a result of the nonlinear evolution. The preservation of multi- 
wavepackets as they evolve shows also that only the nonlinear interactions between small 
neighborhoods of points (k,;, ni) are essential and contribute constructively to the nonlinear 
dynamics, whereas the amplitudes of modes with wavevectors k outside those neighborhoods 
is vanishingly small as /3, g — > 0. The later is quite remarkable since the coupling term 



by rigorous analysis, is based on a fact that the contribution of this term to the solution 
is a time integral involving highly oscillatory functions that becomes vanishingly small as 
(3, g 0. This mechanism is similar to the classical averaging mechanism for systems of 
ordinary differential equations described, for instance, in |11| : the relevance of the averaging 
mechanism for long- wave asymptotics for hyperbolic systems of PDE is well-known, see [30] . 

We would like to relate now the multi-wavepacket preservation property to the linear 
superposition for wavepackets established in [7]. According to that principle if the initial 

state h = ^^hi, with h;, I = 1,...,N being "generic" wavepackets, then the solution 
u{t) = Q (h) (r) to the evolution equation ([TS]) equals with high accuracy to the sum of 
individual solutions u; of N equations with respective initial data h;. Namely, if /3, p > 
satisfy the following relation 



provided proofs of (|^ remain valid as long as 1(1^ holds. Obviously, the bound 

P > CiQ in ([1^ determines when (|44p becomes trivial. This bound is sharp and examples 
below show that when (3 ^ g the remainder D (r) in does not tend to zero when /3 0. 

The both the multi-wavepacket preservation and the linear superposition apply to sums 
of generic wavepackets. It is important no notice though that the multi-wavepacket preser- 
vation holds for any dependence between g and /3 which satisfy that is g{l3) < C(3'^ 
with arbitrary small q whereas the linear superposition holds if £»(/?) < C(3. Thus, the 
bounds ([^^ on (i determine the range of its values for which the both multi-wavepacket 
preservation and the linear superposition hold simultaneously (provided some genericity 
conditions are satisfied). In this range wavepacket preservation provides additional infor- 
mation on behavior of solutions with single wavepacket initial data, namely that the solution 
remains a single wavepacket. Obviously, linear superposition principle does not follow from 
multi-wavepacket invariance. Below we use simple examples and models to discuss different 
ranges of parameters g and fi where wavepacket preservation is valid but the solutions of 
equations exhibit different behavior. 




A qualitative explanation to that, confirmed 



/3, e ^ 0, P > Cig with some Ci > 0, 



(42) 



then for all times < t < r, we have 
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2.1 An exactly solvable model and the effect of large group velocity 

Here we introduce a simple exactly solvable model for our general system ([T]) which makes 
explicit that in the limit p — > nonlinear effects do not vanish, in particular the blow-up time 
does not tend to infinity. This example also shows that on the time scale where r is of order 
1 solutions undergo significant nonlinear evolution. The influence of g on solutions through 
the group velocity in this example can be seen explicitly. The model is the following system 
of two coupled nonlinear first order hyperbolic equations for variables ui [x, t), U2 {x, t) with 
one-dimensional spatial variable x: 

drUi = -—dxUi+ Fi{ui,U2) , (45) 
9 

drU2^ -dxU2+ F2{ui,U2) , Ci ^ C2, | ^^q = ^1 (^) ' ^^1^=0 = ^2 (x) , (46) 

where the initial data hi , /12 in (|46p are of wavepacket form: 

hi {x) = (Px) coski^x, /i2 (a;) = $2 (/3a;) cosfc2*a;, | 7^ |fc2* | . (47) 

We take the nonlinearity to be quadratic and of the following simple form 

Fi{ui,U2) = ul + aiUiU2, F2 (ui, ^2) = W2 + a2MiU2- (48) 

The system (P5|) - (H7|) allows for explicit form of solutions with one- wavepacket initial data, 
describing a wave propagating with a constant speed controlled by the linear part and with 
a shape evolution controlled by the nonlinearity. This simplest case is compared then with 
the case of two-wavepacket initial data, for which explicit solution is not available. 

In the case when h2 — the second equation has trivial solution U2 = and the system 
(|l5l) - (HS)) reduces to a single equation (|l5l) . The solution to this equation has the form 

of a traveling wave vi (^x — ^t, where vi {y, r) is a solution of the ordinary differential 
equation 

drvi = Fi {vi,0) , vi {y, 0) = hi (y) . (49) 
The explicit formula in the case (|49)l yields 



hi(x-^) ^JpL---^))coski,(3(x---f) 

VI (x, r) = ^ ^ \ = } f W / N (50) 

l^rhii^x~^j l-T$i(/?(a;-^jjcosfci,/3(a ' 



CjT_ 

e 



for a time interval < r < tq where tq = sup \hi{y)\ blow-up time. Obviously, the 

blow-up time does not depend on g. Consequently, the wave propagates with the velocity 
^ with its shape evolution being controlled by the nonlinearity. Similarly, when hi = 



the first equation has the trivial solution ui = and the system (|45 |) - (|46|) reduces to a 
single equation (|46p which has a solution in the form of a traveling wave f 2 (^x — ^t, 

propagating with the velocity Observe that for the simple model (|45|) - (|47)l the group 
velocity coincides with the velocity of a traveling wave. 

The above model is not exactly solvable if the both initial conditions hi and /i2 do not 
vanish. But one can still see the way g influences the nonlinear dynamics quite explicitly by 
applying the superposition principle from 6 . Indeed, let us assume that hi and ft.2 are two 
nonzero initial wavepackets. Then the approximate superposition principle is applicable (in 
order to put the system in the framework of [6] we use 4-component extension (I115|) and 
set g = (3^ , x' > 1). According to the principle the exact solution (ui, U2) is approximated 
by (^vi (^x — ^T, ,V2(^x — ^r, r^^ , which is explicitly given by (j50p with the accuracy 

O ( iTrfT ) =0(/3^^^'^) with arbitrary small e if ci / C2. As it as shown in [6J the validity 



.0 . 

of such approximate presentation is due to the large difference ^^^^ of the group velocities 
of two wavepackets. 
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2.2 Dispersive effects and nonlinearity 

Based on an elementary example of the Nonlinear Schrodinger equation (NLS) 

drU = — [70U + i'yidxU + 72^^^] + bi \u'\' u, u — u{x,t) , x (51) 

with the initial data in the form of a wavepacket u\^^q — $ (/3x) e''^*^ we would like to 
explain here why we are interested mostly in the case 

^ > C > 0, (52) 

when the dispersion is not dominant. To make the dependence of w on /3 and g explicit we 
change the variables 

uix) =v{px)e'''''', (3x^z, (53) 

and obtain equation 

drv - [^ov + ^Pl'idzv + -f^f3^dlvi] + b \v\^ V, H.=o = * (^) , (54) 

where "f'l = 7i//3 + 272 fc*. Changing variables once more 

' f P \ P 
V {z, t) = e~~S''^ow I z H — 7'j^T, r j , z H — ■y'^r — y, (55) 

\ Q / Q 

we obtain for the envelope w the following standard NLS equation 

drw = -—j2dyW + b \w\^ w, wL „ = *(?/), < T < r,, (56) 
g 

with initial data independent of the parameters /?, g. The behavior of the solution w to the 
equation (|56p on the time interval < r < is determined by the dispersion parameter 

and evidently linear dispersive effects become significant when ^ is not too large. If 

^ > 00 and /3 ^ the solution tends to zero at every fixed r = tq > 0. Indeed, if we take 

g — , > 2, and make another change of variables r = t[3^ , w = (3^^^ ^'^W , the 



equation (|56p reduces to the following problem with small initial data 

dtW ^ -i-l^dlW + b\W\^W, VKlj^o = /3^'/2-^$(y) . (57) 

For small enough (i the solution W to this problem exists for all t and W {t) {) as t 00 
(see [16 ). In particular, for t = tq/?^^^ we have w (tq) when /? — * 0. 

In the general case, the solution dependence on small /3, g is as follows. The dependence 
on large i in ([51]) is completely described by the change of variables (|55|). yielding a wave 

which (i) moves as a whole with a large group velocity ; (ii) has a slowly evolving shape 
as described by v and w in ([55)) . ([55)l . ([SS]) . 

The above observations show that for small the dispersive effects dominate and control 
the nonlinear ones. Keeping that in mind and being interested in stronger nonlinear effects 
we focus primarily on the case ((52)) . i.e. ^ > C > 0, for which there are two scenarios of 

the nonlinear evolution. In the first scenario, when > 0, the linear dispersion produces 

only a small correction to the solution of the equation drW = b \w\^ w with that nonlinear 
equation governing the nonlinear dynamics of the envelope w for being smaller than the 
blow-up time. In the second scenario, when 0^ ^ g, the equation ([56|l becomes independent 
of /?, g and describes the evolution of the envelope w governed by an interplay between the 
dispersion and the nonlinearity. The case (3^ ^ g can be also characterized as one where 
dispersive effects do occur but they don't dominate nonlinear effects, and, as it is well known, 
the dispersion can exactly balance the nonlinearity yielding solitons. 
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2.3 A coupled modes system 

Here we illustrate statements of the general theorem on the wavepacket preservation and 
the approximate superposition principle by a simple but still nontrivial example. Let us 
consider a system of two coupled NLS type equations for variables Ui {x,t), U2 (x^t) with 
one-dimensional spatial variable x 

= |_7oi + ilii^^x + 72i9^J wi + (^fcii + bi2 \u2\ J ui + C12 \u2\ U2, (58) 

drU2 = [702 +1712^^:1 + 722^x] ""2 + (b21 +622 |U2|^) "2 + C22 l^ll^Ul, (59) 

"i|.=o = '^i (^) = *i (/?^) e"*'", W2|,=o = (2:) = $2 iPx) e•'=-^ (60) 

where 7^ are real and 6y are complex coefficients and the initial data in (j60p are in the form 
of wavepackets with $j (y) being Schwartz functions. Notice that if in the coupled modes 
system (I58l) - (|60p /12 = and C12 = C22 = then it has trivial solution U2 = 0, and reduces 
to a single NLS equation of the form ((5T|) . The dependence of the solution {ui, U2} on the 
large i is captured by the change of variables (|55p . Namely, ui is a wave with a slowly 

varying envelope described by vi which moves with large velocity . The dependence on 
(3 is of the form vi {y, r) = wi {(Jy, t) (see following subsection for details). Similarly we can 
consider the case when hi = for which the first equation has trivial solution ui — 0, so the 
system (|58l) - ([59)) reduces to a single equation (l59l) with the solution represented by a wave 
having large spacial extension proportional to ^ and moving with the large velocity 

2.3.1 The superposition principle 

Let us assume here that hi 7^ 0, /12 ^ 0, C12 ^ 0, C22 ^ and (3 = g'^, < k <1. Applying 
the superposition principle we obtain for generic fc,i, fc*2 the following representation of the 
exact solution 

ui (x, r) = vi {x, r) e'''"^'' + Di, U2 (x, r) = V2 (x, r) e'''-^^ + D2 

where vi (x, t) is a solution of the NLS equation (|58p with hi2 = C12 = 0, with V2 {x, t) being 
a solution to a similar decoupled NLS equation for 622 = C22 = 0, and Di and D2 are small 
terms satisfying 

supo<,<,J|i?i(-,T)||^^ +supo<,<,J|Z?2(-,r)||^3. <Cr'-i-^+C/?, >/ ^ >c-\ (61) 

1 2 I 1 2 

We would like to emphasize here that the coupling terms 612 |u2| ui + C12 \u2\ U2 and 
&21 1^1 1 U2 + C22 \u2\ U2 in the equations |5($|) -/ f5Pj) are not small whereas their ultimate 
contributions to the solutions are small. One can explain/interpret that phenomenon as 
being due to the destructive wave interference and mismatch of group velocities. 

2.3.2 Wavepacket preservation 

Here we assume that hi ^ 0, /12 = 0, C12 ^ 0, C22 ^ and g — j3'^ , < x' < 2. According 
to the wavepacket preservation we have 

Ml {x, t) = vi {x, t) e"'*^'' + Di, U2 (x, r) = Di, 

where vi {x, r) is a solution of (j58p with 612 = , C12 = 0, and Di and D2 are small terms 
satisfying 

suPo<r<r. \\Di {■,t)\\^^ +supo<^<^^ \\D2 {■,t)\\^^ < Cg 
Notice once more (see the above section) an interesting phenomenon: the equation (|59p 
for U2 (x, t) has a coupling term &21 |^i|^ "2 + C22 l^i |^ ui which docs not become small as 
/3, g — + 0, but, remarkably, its ultimate contribution to the solution is small. 
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2.3.3 Limitations of the superposition principle 

Now we provide an example based on the system (j58p - (|60p with Ci2 = C22 = showing that 
the above estimate (|61l) in the superposition principle is sharp in the sense that -i-e 
cannot be replaced by -i+<^ with V > 1. We set here k' = \ and g = f3. After the 
change of variables ([55]) for mi,U2 followed by yet another change of variables /3x = z, 

Vl = e '^^wi, W2 = e '^t~W2 we obtain from ([55]) - ([50)1 the following system: 
drWi = -i [ij'iidzWi + Pj2i9'^u}i] + (bn \wif + 612 |w2|^) wi, 

l9rW2 = -i [i7i2f^2W^2 + /3722'9z^"2] + (^^21 + ^22 1^21^) W2, 

«^llr=0=*l (^)' ^^2U^o = *2 (z). 

This system has a regular dependence on /3 as /3 — > with the solution converging in L°° 
to the solution of the system with /3 = 0. If we set now in the last system 612 = 621 = 
it turns into a system of two decoupled equations. Notice then that the difference between 
the solutions of the decoupled system and the original one does not tend to zero as /3 ^ 0, 
implying that the superposition principle does not hold when q = j3. 

2.4 Wavepacket interaction system with a universal nonlinearity 

We will prove in the following sections that the dynamics of a multi-wavepacket with a 
universally resonance invariant nfc-spectrum for a general system can be approximated with 
the accuracy 0{q) by substituting the nonlinearity with a properly constructed universal 
or weakly universal one. Here we provide an example of a system, called wavepacket inter- 
action system^ with a universal nonlinearity and show that its dynamics preserves simple 
wavepackets as in (jl2l) . It is shown later that universal nonlinearities are related to univer- 
sally invariant multi- wavepackets in the sense of Definition [T8l 

Wavepacket interaction system with universal nonlinearity has the form similar to NLS, 
namely 

drUj^C = - [-iC7oj + 7i J • VrM,,c - iCVr • 72j Vr%,c] + ^j,C (") > I- e M'*, (62) 
Q 

u = . . . , UN+, UN-) , j = 1, . . . , iV, C = ±, (63) 

".,cl.=o = h^.C, h,.^ (r) - (/3r) e'^^-^, (64) 

where for every j coefficient 7q j G R, j £ M.'^ is a vector, 72 j is a symmetric dxd matrix, 
lij • Vr is a first order scalar differential operator, Vr • 72 jVr is the second order scalar 
differential operator, and the universal polynomial nonlinearities have the following form: 

FjX (u) = ^^J^C 11;=! (UL + UL-)"' UjX, (65) 

where P = (t^i, • . • , i^n) , j = 1, . . . , A^, C = ±- 

Remark 8 Notice that if we set hj,^ = h*_^, bp_j^-^- — j - ^ ^1^,3 '^"■'^ ~ - ~ 
then ui,+ui^- = \ui^^\^ and Fj^^ (u) turns into 

Fj {ui,...,un) X^^liX^ipi^^^^jIli^i 1"'!^''' ""j"' (^^) 
and equations of i62\) with ( = + turn into 

drU-j = i [-i7oj + llj ■ VrUj - iVr • 72jVrWj] + Fj{ui,..., un) , (67) 
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Obviously, a solution of |6'7[ ) defines a solution Uj^+ = Uj, Uj_^ ~ u* of the system 
In the simplest case N — 1, d — 1 ^67^ takes the form of classical NLS: drU = -^dxU ~ 

i^dlu + b\u\^u. 

Note that the universal nonlinearity Fj,(^ has a characteristic property 

(68) 

holding for arbitrary set values (j)^. We also consider more general nonlinearities F for 
which ()68p holds for a fixed set of fi'equencies 0; — ujm (k*;), and call them weakly universal. 
We introduce now the averaging operator At acting on polynomial functions F : (C^)^ — > 
(C^)"" by 

(^Tn-,C=K>^)^,,= (69) 

1 F 

- / e-";*.*F_,- c (e''^i*ui,+ , e-'^i*?/i,_, . . . , e'^"*ujv,+, b-'^-'un.-) dt, 
^ Jo 

where cf) — {(j>i , . . . , (/'at) • The operator Aj, ^ depends on the frequency vector (j) — {(j>i , . . . , (f)]^) 
If is a universal polynomial nonlinearity, then (^Aj, = Fj^i^ for any choice of fre- 

quencies (jji, . . . , (j)fq. Note that averaging 



G'av,,,c (^) = lim (a^ -^G) {u) (70) 

is defined for any polynomial nonlinearity G : (C^)^ (C^)^. If is generic, then 
Gavjx i''^) always a universal nonlinearity. In a general case Ga^v.j,^ for given frequencies 
4> one obtains a weakly universal nonlinearity which might be not universal. 

Systems with universal nonlinearities have interesting properties which we describe in 
the following Proposition and remark. 

Proposition 9 Let g ~ P and ^2 j = 0- Then evolution governed by the first order system 
with a universal nonlinearity I16S\) preserves simple wavepackets as defined by il2\) . 



Proof. Let u(t) be a solution of (p^ for < r < r*. Using the property 
change variables 



and obtain from ([62 



drVjx = ^7i, ■ ^rVj.c + F,,^ (v) , vj,^l^^ = ^ (/3r) . (72) 
Changing variables 

Vj,c (r, t) = Wjx (/3r, t) , /3r = z, (73) 
we obtain from (l72|) that wj is a solution of the following system of differential equations 

which does not depend on f3. Then using (|73p and (I7ip we observe that every component 
ui of the solution to (j62p has the form of a simple wavepacket for every r G [0,t»], with an 
envelope lij (r). ■ 
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Remark 10 Equations i62\) with universal nonlinearities allow special solutions in the form 

of Uj^i^ ~ e^^*^''^e~^~'^ Vj^c^ (r) where Vj^(_ (r) do not depend on r. // the initial data in {IS^ 
are constants, ^j^t^ (/3r) = '^j^Q (0), then \7i3^ turns into a system of ODE. This implies that 
every linear suhspace of pure modal functions with the basis Uj-e**^'', Uj__e~'''*j"", j = 1, .., N 
is invariant with respect to nonlinear equations i6S\) . Another class of special solutions of 
\6'2\) are time-harmonic solutions of the form Uj,^ (r, t) = e~'^'^^'^Uj_^ (r) where Vj^^ solve a 
nonlinear eigenvalue problem; for universal nonlinearities LUj can be considered as unknown 
nonlinear eigenvalue. Existence of such special solutions is a special property of universal 
and weakly universal nonlinearities. It is remarkable that original nonlinear equations might 
not have time harmonic solutions whereas equations with canonical nonlinearities which 
approximate evolution of wavepackets (see Theorem^ admit such solutions. 

2.5 Invariance of excited modes for finite-dimensional ODE's 

Here we discuss the resonance invariance conditions imposed in Theorem [5] in a simpler case 
of finite-dimensional ODE's. In this case one can also see the rise of universal nonlinearities 
in the process of time averaging. As we already discussed in the introduction, a PDE system 
^ when restricted to constant functions turns into the following system of ODE's 

drV = --L"U + F (U) , U (r)l^^o = h, h e C^^, U e C^-', (75) 

where F (U) is a polynomial, U — {Ui^+, Ui^-, . . . , Uj^+, Uj-) G C^"'. We assume that the 
eigenvalues 0Jn,c, (0) — uj'^ ^oi the Hermitian matrix L*' = L (k)|jj.^Q are distinct 7^ + 
for j ^ i and the symmetry conditions ^ take the form lo^ = — ^. We also assume 
that the eigenvectors of L*^ coincide with the coordinate orts in C^'^. The following limit 
case of Theorem [5] with = shows that solutions to this system have the property to 
preserve the set of initially excited modes. 

Theorem 11 Let the initial data h = . . . ,hj,^,hj,^) G C^"^ in |75[ j have non- 

zero components hj^c^ only for a subset B of indices j € {1, . . . , J}, and let B' — {1, . . . , J} \ 
B be its complementary set. Assume that B is resonance invariant in the sense that the 
resonance equation 

""n'X - T^Zi """..C"' = ^' ^'^^^^ ^ ^' ^^'^ + (76) 

does not have solutions if n' E B' ( compare with Definition [I^ in the special case when all 
k,; = Oj. Then under the nonlinear evolution of |75[ ) modes with indices n' G B' remain 
essentially unexcited in the following sense 

sup \Un' (t)| < Cg for all n' e B' . (77) 

0<T<T. 

Note that F (U) provides a nonlinear coupling between modes U^. ^(j) with nj G B 
and Un'x with n' G B', but the resulting interaction is not 0(1) on a fixed time interval 
[0,T,] as one might expect, but rather of order O (v) as (|77p shows. One way to prove 
Theorem [TT] is to follow the proofs of Theorems [35] and [37] with obvious modifications and 
simplifications. In particular, instead of (jlSp one has to consider the following system with 
oscillatory coefficients 

a^u = ef ^"f (e^^"u) , u(T)|^^o = h. (78) 

Alternatively, Theorem[Tl]can be derived directly from the classical time averaging principle. 
Indeed, the time averaging of (|78p yields the following averaged system 

9.v = Fa.(v), v(r)|^^o = h. 
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where Fav is defined as in ([M]) . fTO)) with the frequencies = From the Krylov- 

Bogolyubov averaging theorem (see [11], [37]) one obtains 

|v(r) - u(t)| < Cp, < T < T,. 

A straightforward examination shows that if B is resonance invariant and j G B' then the 
polynomial components Favj,c ("^) factorize into Fav.j,c C'^) ~ X^j'es' c' ^av j' c' 
implying ([77)1 since Uj^^ (0) = for j & B' . 

A stronger universal resonance invariance condition in Definition [TH] also takes a simpler 
form in the ODE case. Indeed, let us collect the terms in ((76|) at different a;°_|_ as in (|101[) . 
namely 

^n'x - 51,^1 ^^^1+^ '^h'^re 6^ are integers, (79) 

Similarly to Definition [15] we call B universally resonance invariant if every solution to the 
resonance equation (1751) must have n' G B and every coefficient 6i in ([75)1 for the solution 
is zero, i.e. 5i = 0, i = 1, . . . , J. Obviously, if all oj" ^ are rationally independent then it is 
universally resonance invariant. 

Now let us look how universal nonlinearities arise under time averaging. Observe that if 
the entire set {1, . . . , J} is universally resonance invariant and (v) are arbitrary poly- 
nomials, then the polynomials F^vjx i^) obtained by discarding the "resonant" terms 
in e 8 Flee uj yielding universal form (|d5|) . (|d6|) . For example, if F is an arbitrary 

cubic nonlinearity in then the time averaging yields NLS-like nonlinearity Fav with 
components 

When B is resonance invariant but not universally resonance invariant the averaging 
produces a weakly universal nonlinearity. A nonlinearity which is weakly universal but 
not universal may include additional terms, for example the cubic nonlinearity in classical 
four-wave interaction system where it is assumed that a;^ - + ^3 + + ^4 + — '-^1 + (see 
[46j p. 201) in the equation for ui^+ in addition to NLS-like terms involves the product 

3 Conditions and definitions 

In this section we formulate and discuss definitions and conditions under which we study 
the nonlinear evolutionary system (|T]) through its modal, Fourier form Most of the 
conditions and definitions are naturally formulated for the modal form ([3]), and this is one 
of the reasons we use it as the basic form. 

3.1 Linear part 

The basic properties of the linear part L (k) of the system ([3]), which is a 2 J x 2 J Hermitian 
matrix with eigenvalues LOn.c, (k), has been already discussed in the Introduction. To account 
for all needed properties of L (k) we define the singular set of points k. 

Definition 12 (band-crossing points) We call ko a band-crossing point for L (k) if 
ujn+i,c, (ko) = '^n,c, (ko) for some n,C, or h (k) is not continuous at kg or i/ wi.± (kp) = 0, 
we denote the set of such points by tibc- 

In the next Condition we collect all constraints imposed on the linear operator L (k) . 
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Condition 13 (linear part) The linear part L (k) of the system (0) is a2Jx2J Hermitian 
matrix with eigenvalues uJnX (^) ^'^'^ corresponding eigenvectors gn.c 0^) satisfying for k ^ 
CTftc the basic relations In addition to that we assume: 

(i) the set of hand-crossing points a be is a closed, nowhere dense set in R'' and has zero 
Lebesgue measure; 

(a) the entries of the Hermitian matrix L (k) are infinitely differentiable in k for all k ^ 
CTfjc that readily implies via the spectral theory, f35}j . infinite differentiability of all 
eigenvalues a;„ (k) in k for all k ^ a; 

(Hi) L (k) satisfies polynomial bound 

||L (k) II < C (1 + \kf) , keW^, for some OQ and p>0. (80) 



Remark 14 (dispersion relations symmetry) The symmetry condition ^ on the dis- 
persion relations naturally arise in many physical problems, for example Maxwell equations 
in periodic media, see l]]l-J^, J^, or when L (k) originates from a Hamiltonian. We would 
like to stress that this symmetry conditions are not imposed to simplify studies hut rather 
to take into account fundamental symmetries of physical media. In fact, the opposite case 
when ( ^ is assumed not to hold is much simpler. The symmetry creates resonant nonlin- 
ear interactions, which makes studies more intricate. Interestingly, many problems without 
symmetries can be put into the framework with symmetry by an extension of the relevant 
system (see Section 4). 

Remark 15 (band-crossing points) Band-crossing points are discussed in more details 
in JT, Section 5.4], 12, Sections 4-1, 4-2]. In particular, generically the set abc of band- 
crossing point is a manifold of the dimension d — 2. Notice, that there is an natural ambi- 
guity in the definition of the normalized eigenvectors gnX (M) L (k) which is defined up 
to a complex number ^ with \^\ = 1. This ambiguity may not allow an eigenvector gnX 0^) 
which can be a locally smooth function in h to be a uniquely defined continuous function in 
k globally for all k ^ a^c because of a possibility of branching. But, importantly, the orthog- 
onal projector Tln.(^ (k) on g„ (k) as defined by ill]) is uniquely defined and, consequently, 
infinitely differentiable in k via the spectral theory, fSSf . for all k ^ abc- Since we consider 
U (k) as an element of the space and abc is of zero Lebesgue measure considering k ^ abc 
is sufficient for us. 

We introduce for vectors u G C'^'' their expansion with respect to the orthonormal 
basis{g„,<; (k)}: 

J J 
" (k) = E E ""'C (k) gn,c (k) = E E ""'C (k) , Cin,c (k) = n„,c (k) u (k) (81) 

n=l(=± n=lC=± 

and we refer to it as the modal decomposition of u (k) and to iinX (k) as the modal coefficients 
of u (k) . Evidently 

'y^ '' ^ ^ n„ (k) — I2J, where I2J is the 2 J x 2 J identity matrix. (82) 

Notice that in view of the polynomial bound [80]) we can define the action of the operator 
L (— iVr) on any Schwartz function Y (r) by the formula 

L (— iVr) Y (k) = L (k) Y (k) . where the order of L does not exceed p. (83) 

In a special case when all the entries of L (k) are polynomials ([55)1 turns into the action of 
the differential operator with constant coefficients of order not exceeding p. 
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3.2 Nonlinear part 

The nonlinear term _F in ([S]) is assumed to be a general functional polynomial of the form 

F ( tj J — ^(™) j U™ J ^ where F^™) is m-homogeneous polylinear operator, 

(84) 

rSOTf}- (85) 

The integer mp (|85p is called the degree of the functional polynomial F. For instance, if 
dyip = {2} or OJli? = {3} the polynomial F is respectively homogeneous quadratic or cubic. 
Every m-linear operator in is assumed to be of the form of a convolution 



dJlp ~ {'Til, . . . ,mp} C {2,3, . . .} is a finite set, and mp = max{r 



F 



^(m) ^ 



Ui, 



,Um) (k,r) 
where I 



^ X^"^ (k, k) Ui (k') . . . U,n (k('") (k, k)) (86) 



dk' . . . dk("-l) 



k(™) (k, A?) = k - k' - . . . - k(™-i), fc = (k', . . . , k(™) 



(87) 



indicating that the nonlinear operator F^"^^ (Ui, . . . ,!!„) is translation invariant (it may 
be local or non-local). The quantities ((86|) are called susceptibilities. For numerous 

examples of nonlinearities of the form similar to (|84p. ([86|) see [I]- [2] and references therein. 
In what follows the nonlinear term _F in ([3|) will satisfy the following conditions. 



Condition 16 (nonlinearity) The nonlinearity F yJj is assumed to he of the form {84^- 

Ii80\) . The susceptibility x''™'' (k, k', . . . ,k(™)) is infinitely differentiate for all k and k*^-'^ 
which are not band-crossing points, and is hounded, namely 



(m) 



(27r 



-(m-l)d 



sup 

k,k',...,k<")GK''\crbc 



Am) 



<Cy, me Mp, 



where the norm 
defined by 



X 



(k,fc) 



(k,k',...,kM 

l^k, fc^ of m-linear tensor x'^™'^ ■ (C^'^)™ (C^-^)™ for fixed \i,k is 
(k,^) (xi, . . . , x„i) , where |x| is the Euclidean norm. (89) 



sup 

Xj|<l 



Am) 



When x*-™-* (^k, fc j depend on small or, more generally, on g'^, q > 0, we simi- 
larly have x*-™-* k, g''^ . Many results of this paper extend to this case, in particular 
(k, k, gi^ - x'"^ (k, k, < C'^g'' for < 1 then conditions of Corollary [35] are 



if 



(rn) 



fulfilled. 

Note that since the tensors x*-™'' (^k, kj are bounded, the dependence on ^k, kj cannot be 

polynomial, therefore the original equation ^ does not include spatial derivatives but rather 
includes bounded " pseudodifferential" operators. Note that this type of susceptibilities with 
spatial dispersion is common in nonlinear optics, see [15j . |41j . [55] . 

3.3 Resonance invariant nA;-spectrum 

In this section, relying on given dispersion relations cj„ (k) > 0, n G {1, J}, we consider 
resonance properties of nfc-spectra S and the corresponding fc-spectra Ks as defined in 
Definition [21 i.e. 

S = {{ni,Ki), ; = l,...,iV}cI] = {l,...,J}xR^ Ks ={k„^, z = 1, |i^s|} • (90) 
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We precede the formal description of the resonance invariance (see Definition I18p with the 
foUowing guiding physical picture. Initially at r = the wave is a multi-wavepacket com- 
posed of modes from a small vicinity of the nfc-spectrum S. As the wave evolves according to 
([5]) the polynomial nonlinearity inevitably involves a larger set of modes [5']^,^,^ 3 S, but not 
all modes in [S'Jq^^. are "equal" in developing significant amplitudes. The qualitative picture 
is that whenever certain interaction phase function (see (|134p below) is not zero, the fast 
time oscillations weaken effective nonlinear mode interaction and the energy transfer from 
the original modes in S to relevant modes from [S'l^y^., keeping their magnitudes vanishingly 
small as P, Q ^ 0. There is a smaller set of modes [Sf^^^ which can interact with modes 
from S rather effectively and develop significant amplitudes. Now, 

if [Sy^^^ C S then S is called resonance invariant. (91) 

In simpler situations the resonance invariance conditions turns into the well-known in non- 
linear optics phase and frequency matching conditions. For instance, if S contains (no, 'k-*io) 
and the dispersion relations allow for the second harmonic generation in another band ni 
so that 2uJno 0^*io) — (2k*/(j), then for S to be resonance invariant it must contain 
(ni,2k,;(,) too. 

Let us turn now to the rigorous constructions. First we introduce necessary notations. 
Let m > 2 be an integer, I — {li, .., l„i), Ij £ {1, N} be an integer vector from {1, A^}™ 

and C = (C^^\,-,C^™^), C^^^ e {+1,-1} be a binary vector from {+1,-1}". Note 
that a pair ^C, tj naturally labels a sample string of the length m composed of elements 
ni^, from the set {+1, —1} x S. Let us introduce the sets 

A = {(C, l):le {1, A^} , C e {+1, -1}} , (92) 
A"' = |a = (Ai, Xm) , Aj e A, j = 1, m| . 

There is a natural one-to-one correspondence between A™ and { — 1,1}™ x {1,..., N}"^ and 
we will write, exploiting this correspondence 

A= ((C',/i),...,(c'™\;™)) ^{It), de{-lAr, I e {!,..., Nr for A G A™. (93) 
Let us introduce the following linear combination 

Xrn (a) = >im (C, t) = , C^'^k,;^. with C''''^ G { + 1, -1} , (94) 

and let ['5']^o^t be the set of all its values as k*;^ G Ks, A G A™, namely 

We call [S]j^ output k-spectrum of Ks- Everywhere in this paper we consider nfc-spectra 
S which satisfy the following condition 

[S]K,ontr\''bc^^- (96) 

We also define the output nk-spectrum of S by 

[^lout = {('^^k) G {1, J} X R"^ : n G {1, J} , k G . (97) 

We introduce the following functions 

^i,m{xj (kjj ^^^_^C''^^^h , fc* = (k*i, ...,k»|Ks|) , where k,i. e Ks, (98) 
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SI (^C,n, )Cj ^k,,, fc,^ = -C^n (k„) + rii,,„ (ji^ ^fc*^ , (99) 

where ( = ±1, m & Tip as in ([51]) . We introduce these functions to apply later to phase 
functions (fTM)) . 

Now we introduce the resonance equation 



n 



(C, n, a) (c>^™ (a) , £) = 0, fe {1, Nr,Ce {-1, ir , (100) 



denoting by P {S) the set of its solutions ym, n, Xj . Such a solution is called S -internal 
if ^ 

(^n, C>^„J (^A^^ e S*, that is n = nj^, (^Cm (xj = Kio, h e {1, -^N] , 



and we denote the corresponding Iq — I \ \j . We also denote by Pint {S) C P (S) the set of 
all 5'- internal solutions to (jlOOp . 

Now we consider the simplest solutions to (jlOOp which play an important role. Keeping 
in mind that the string I can contain several copies of a single value I, we can recast the sum 
in IMl) as follows: 



(a) = n,,^ (C, l) = ("^^'^ ' "'^^^^ ^' = { ^'''"o*'^ ^^'^ !f (!) i ' 

(101) 

r-i (0 = {j -.1,^1, 1 < i < m} , r= (/i, ...,;„), 1 < ? < A^. 

Let us call a solution |^m, C, A^ G P (S*) of pOOp universal if it has the following properties: 

(i) only a single coefficient out of all 5i in (|lGip is nonzero, namely for some /q we have 
= ±1 and (5/ = for I ^ Iq; (ii) n = n/^ and ^ = (5/^. yl justification for calling such a 
solution universal comes from the fact that if it is a solution for one A:* it is a solution for 
any other fc* € M''. We denote the set of universal solutions to (|100p by Puniv {S), and note 

that a universal solution is a S'-internal solution with / ^A^ — Iq implying 

Puniv (S) C Pint (S) . (102) 



Indeed, observe that for Si as in (|101l 

(a) = X,n (C, t) = = EL 

implying >im (j^^ = '5/nk*/o and C^m(>^^ = (^/^k^/^ = k»/j,. Then equation (jlOOp is 
obviously satisfied and (n,(xm i^)) — ("-/oik,/,,) £ 5. 



Example 17 (Universal solutions) Suppose there is just a single band, i.e. J = 1, 
a symmetric dispersion relation uji (— k) = cui (k), a cubic nonlinearity F with ^Xftp — 
{3}. First let us take the simplest nk-spectrum Si = {(l,k*)}, that is N — I. Then 

^1.3 ^fc*^ = SiLOi (k*) and Km = (5ik* where we use notation \101\) . The uni- 
versal solution set has the form Puniv {Si) = | ^3, 1, : A G A^, = ±| where A+ con- 
sists of vectors (Ai, A2, A3) of the form ((-, 1) , (+, 1) , (+, 1)), ((+, 1) , (-, 1) , (+, 1)) and 
((+, 1) , (+, 1) , (— , 1)). Obviously, Puniv{Si) — Pint{Si). In the next example we take the 
nk-spectrum S = {(l,k:„) , (1,— k^^)}, that is N = 2 and k*i = k*,k,t2 = — k:„. This example 

is typical for two counterpropagating waves. Then ^1^3 (^Xj = C'"'''^ij — 
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{Si + S2)u!i (k,) and k^. ^Aj = X^jli C*'"'^k,;^. = Sik^i +82^*2 = (Si — 62) where we use 

notation UOl]) . The universal solution set has the form Puniv [S] — | ^3, 1, : A G A,^, C = ± | 

where A+ consists of vectors (Ai, A2, A3) of the form ((+, 1) , (— , 1) , (+, 1)), ((+, 1) , (— , 1) , (+, 2)), 
((+, 2) , (— , 2) , (+, 1)), ((+, 2) , (— , 2) , (+, 2)), and vectors obtained from the listed ones by 
permutations of coordinates Ai, A2, A3. The solutions from Pint {S) have to satisfy \5i — 52\ = 
1 and |(5i + (52| = 1 which is possible only if 5i52 = 0. Since C, ^ 5i + 82 we have 
C^Km {t^i = [5l - Sf) k, and C,Xm, [xj = k,i if \Si\ = 1 or C^^m [pj = k,2 if |^2| = 1- 
Hence Pint {S) = Puniv {S) in this case. Note that if we set S2 — {(1, — k,)} then 5 = 5*1 U 6*2 
but Pint{S) is larger than Pint (Si) U Pint(.S2)- This can be interpreted as follows. When 
only modes from Si are excited, the modes from S2 remain non-excited. But when the both 
Si and 5*2 are excited, there is a resonance effect of Si onto S2, represented, for example, 
by X = ((+, 1) , ( — , 1) , (+, 2)), which involves the mode C^m (j^ = k*2- 

Now we are ready to define resonance invariant spectra. First, we introduce a subset 
[5'] out of [S]out by the formula 

[S]'Z = {(", k„) e [Sl^, : k,, = C><m (a) , m e Tip, where (104) 
(^m, n, Xj is a solution of pOOp j , 
calling it resonant output spectrum of S, and then we define 

resonance selection operation TZ {S} = 5 U [SY°^^ . (105) 

Definition 18 (resonance invariant rifc-spectrum) The nk-spectrum S is called reso- 
nance invariant if TZ{S) = S or, equivalently, [-SJ^^'j C S. The nk-spectrum S is called 
universally resonance invariant ifTZ (S) =^ S and Puniv (S) = Pint (S). 

It is worth noticing that even when a nfc-spectruni is not resonance invariant often it 
can be easily extended to a resonance invariant one. Namely, if TZ^ (S) H (Tbc = for all j 
then the set 

n°° (S) = I f (5) C S = {1, J} X R'' 
^j=i 

is resonance invariant. In addition to that, TZ°° (S) is always at most countable. Usually 
it is finite i.e. TZ°° (S) — TV [S] for a finite p, see examples below and we also show below 
that 7^°° {S) = S for generic Kg. 

Example 19 (resonance invariant nfc-spectra for quadratic nonlinearity) Suppose 
there is a single band, i.e. J ~ 1, with a symmetric dispersion relation, and a quadratic non- 
linearity F, that is = {2}. Let us assume that k, ^ 0, k,,2k,,0 are not band-crossing 
points and look at two examples. First, suppose that 2lo\ (k*) 7^ lo\ (2k*) (no second har- 
monic generation) and loi (0) ^ 0. Let us set the nk-spectrum to be the set Si — {(l,k,)}, 
then Si is resonance invariant. Indeed, Ksi = {k,}, [Si\j^ — {0, 2k,, — 2k,}, [^ij^^^ — 
{(1, 0) , (1, 2k,) , (1, — 2k,)} and an elementary examination shows that ['S'l]^^''^ — C Si 
implying TZ{Si) = Si. For the second example let us assume uji (0) ^ and 2uji (k,) — 
wi (2k,), that is the second harmonic generation is allowed. Here [Si^^^^ — {(1,2k,)} 
and TZ{Si) = {(1, k,) , (1, 2k,)} implying TZ{Si) ^ Si and, hence, Si is not resonance 
invariant. Suppose now that 4k,, 3k, ^ abc and wi (0) 0, wi (4k,) ^ 2^1 (2k,), 
oji (3k,) ^ uJi (k,) +LUi (2k,) and let us set S2 — {(l,k,) , (1,2k,)}. An elementary exami- 
nation shows that S2 is resonance invariant. Note that S2 can be obtained by iterating the 
resonance selection operator, namely S2 = TZ (TZ (5*1)). Note also that Puniv {S2) 7^ Pint {S2). 
Notice that oJi (0) — is a special case since k = is a band-crossing point, and it requires 
a special treatment. 
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Example 20 (resonance invariant nfc-spectra for cubic nonlinearity) Let us consider 
one-band case with symmetric dispersion relation and a cubic nonlinearity that is dJtp — {3}. 
First we take Si = {(l,k*)}, we assume that k*,3k* are not band-crossing points, implying 

['5'i]/f,ont = {k*, -k,, 3k*, -3k*}. We have fii^a ^A^ ^fc*) = Z]j=i C^^''^i (k*) = (Ji^i (k*) 
and x-m ^A^ = (5ik, where we use notation \101]) . Si takes values 1, —1, 3, —3 . Ifiuii (k») ^ 

LOi (3k,) then U00\) has a solution only if \5i\ — 1 and Si — C, hence C^rn — k, and 

every solution is internal. Therefore, [SiY^^^ — and TZ{Si) = 5*1. Now consider the 
case associated with the third harmonic generation, namely 3uJi (k,) — toi (3k,) and as- 
sume that LUi (3k,) + 2(jJi (k,) ^ cui (5k,), 3u!i (3k,) ^ uji (9k,), 2wi (3k,) + uji (k,) ^ 
ui (7k,), 2uji (3k,) — LUi (k,) ^ lui (5k,). An elementary examination shows that the set 
S4 = {(1, 3k,) , (1, k,) , (1, — k,) (1, — 3k,)} satisfies Tl{S4) = S4. Consequently, a multi- 
wavepacket having 84 as its resonance invariant nk- spectrum involves the third harmonic 
generation and, according to Theorem\^ it is preserved under nonlinear evolution. 

The above examples indicate that in simple cases the conditions on k, which can make 
S non-invariant with respect to TZ have a form of several algebraic equations, therefore, for 
almost all k, such spectra 5* are resonance invariant. The examples also show that if we 
fix S and dispersion relations then we can include S in larger spectrum S' — TV (S) using 
repeated application of the operation TZ to 5", and often the resulting extended nfc-spectrum 
5" is resonance invariant. We show in the following section that nfc-spectrum S with generic 
Ks is universally resonance invariant. 

Note that the concept of resonance invariant nfc-spectrum gives a mathematical de- 
scription of such fundamental concepts of nonlinear optics as phase matching, frequency 
matching, four wave interaction in cubic media and three wave interaction in quadratic me- 
dia. If a multi-wavepacket has a resonance invariant spectrum, all these phenomena may 
take place in the internal dynamics of the multi-wavepacket, but do not lead to resonant 
interactions with continuum of all remaining modes. 

3.4 Genericity of the nfc-spectrum invariance condition 

In simpler situations, when the number of bands J and wavepackets N are not too large, the 
resonance invariance of nk- spectrum can be easily verified as above in Examples 1 1 9[ [20l but 
what one can say if J or are large, or if the dispersion relations are not explicitly given? 
We show below that in properly defined non-degenerate cases a small variation of Kg makes 
S universally resonance invariant, i.e. the resonance invariance is a generic phenomenon.. 
Assume that the dispersion relations a;„ (k) > 0, n £ {1, J} are given. Observe then 

that rim ^Ci '^j ^ ~ ^rn (^Ci ^ k,|^g|) defined by (|99p is a continuous function of 

k,( ^ CTfjc for every m, C, n, A. 

Definition 21 (cj-degenerate dispersion relations) We call dispersion relations ujn (k), 
n = 1, J, Lo- degenerate if there exists such a point k, G R'' \iTbc that for all k in a neigh- 
borhood o/k, at least one of the following four conditions holds: (i) the relations are linearly 
dependent, namely X]n=o ^"^n (k) ~ '^O: where all Cn are integers, one of which is nonzero, 
and the Cq is a constant; (ii) at least one of aj„ (k) is a linear function; ( Hi ) at least one of 
LOn (k) satisfies equation Cw„ (k) — ujn (Ck) with some n and integer C ±1; (iv) at least 
one of uJn (k) satisfies equation ujn (k) = uJn' (— k) where n' ^ n. 

Note that fulfillment of any of the four conditions in Definition [21] makes impossible 
turning some non resonance invariant sets into resonance invariant ones by a variation of 
k,i. For instance, if 3JIf = {2} as in Example fT9l and 2loi (k) = lui (2k) for all k in an open 
set G then the set {(l,k,)} with k, G G cannot be made resonance invariant by a small 
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variation of k:». Below we show that if dispersion relations are not w-degenerate, then a 
small variation of k»; turns non resonance invariant sets into resonance invariant. 

Theorem 22 If Vim {^C,,ni^,\^ ^k'^j^, k'^|^^,|^ = on a cylinder G in (M"* \ dhc) which 

is a product of small halls Gi C (K"^ \ (Tftc) then either ^m, C, no, G Puniv {S) or dispersive 
relations cOn (k) are LO-degenerate as in Definition \21\ 

Proof. Collecting similar terms in (jlOOp we obtain the following equation for k^ from 

1X1' f (kO = Ct^no ( V'. ff^^kj where (5 (5- are integers. (106) 

Comparing (|106p with (jlOip we see that may be non-zero only if (n, k^) G 5*, that is 
(n,ki) = (n;,k;) with I G {1, A^}, where I = I {i,n) is uniquely determined and (5^„ = Si 
with (5; as in (|10ip . If there are two nonzero coefhcients Si in p06p we use an elementary 
Proposition [24] below noticing that we are in the case (ii) of Definition [2TJ If we do not 
have two nonzero S'^ then either all S'l = or only one S[ = S[^_^ 7^ 0. If all S[ = then the 
right-hand side of (|106p turns into ojng (0) and, Gi C (M.'^ \ ate) , t^ng (0) 7^ 0. Hence, for 
every i the sum X]n=i Kn'^n (ki) is constant, one of (5^„ is non-zero and we are in the case 
(i) of Definition [2TJ If only one S'i^Q with i = io we have 

, E'"", ' (^») = ^'^"o (-^lok'o) for aU k, e G„ k,„ e G,„, (107) 

implying linear dependence of the dispersion relations, namely 

Si ijJn (ki) = Gi, i 7^ ioi where Gi are constant. 

n— 1 

The above equations would not imply linear dependence as in case (i) of Definition [21] only 
if 

C=0, z^zo, n=l,...,J, (108) 
and in this case the equality (|107p takes the form 

E„=i ^'^on^^ C'-'o) = C^^»o ('^^o'^'o) for all k,„ e GiQ. (109) 

Note that in this case we deduce from ([94| and (|98)) that Ylin=i ^'ion — ^'io- \^'io \ 7^ 1 '^e 
are in the case (iii) of Definition I21i whereas if |(Sijj| =1 and n ^ uq wc are in the case (iv) 
of Definition [2TI If = 1 and n = no (|109p turns into Si^ujng (ki^) — Qu^no (<^io^io)- Since 
Ulna > it imphes S[^^ = ( and Uno O^to) = ^no {CKo^io)- Hence, in this case (m, C, ^o, e 
-Puniv (<5'), and since all possibilities are exhausted the proof is complete. ■ 

Theorem 23 (genericity of resonance invariance) Assume that dispersive relations ujn (k) 
are continuous and not Lu-degenerate as in Definition \21\ Let JCrinv be a set of points 
(k,ti, ...,k^|^g|) such that there exists a universally resonance invariant nk-spectrum S for 
which its k-spectrum Ks = {k*i, k^|/fg| } . Then JCrinv is open and everywhere dense set 

Proof. The fact that /Crinv is open follows from the Definition [TS] and the continuity in 
k of the dispersion relations cj„ (k). Let G be a small open ball such that its closure G C 

(E** \ (Tfec) ' ' • It suffices to prove that G fl A^rinv contains at least one point (k* 1 , . . . , k^ | | ) . 
For a given finite set VJlp let us consider all possible 



(m,C,no,A)eU x{-l,l}x{l,...,J}x 
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which are not universal solutions to (jlOOp . and for a given yn, no, Aj let Go yn, ^, uq, Xj 
be a set of solutions (ki, ...,k|if,|) to (jlOOp in G, and notice that it is a closed set. Let now 
Go (S) C G be the union of the sets Gq (^m, (, uq, over all ^m, (, uq, A^ £ P (<S')\Puniv (S) 
and let us show that Gq {S) ^ G. Indeed, suppose that Gq {S) = G and hence G is a finite 
union of closed sets. According to the Baire's theorem one of the sets Go (^m, C, uq, Xj with 

(yn, C, no, Xj ^ P {S)\Puniv (S) must have a nonempty interior. Then, according to Theorem 

[22l the dispersion relations ujn (k) are w-degenerate as in Definition [21] contradicting the 
conditions of the theorem. Hence, there is always a point (k*i, k^j^^i) G P (5) \Puniv {S) 
that completes the proof. ■ 

The proof of the next statement is elementary and we skip it. 

Proposition 24 Let fi (k), /2 (k), (k) be real-valued and continuous functions respec- 
tively in neighborhoods of k*i, k*2, k*i + k*2 inMf^. Assume that the following equation 

/l (ki) + /2 (k2) = /3 (<5lki + <52k2) + Go 

holds in these neighborhoods where Co,Si,S2 are constants and S162 ^ 0. Then all three 
functions fi (k), /2 (k), f^ (k) are linear in neighborhoods o/k*i, ]i^2, k*i+k*2 respectively. 



4 Reduction to a standard framework 

Many well known nonlinear evolutionary equations and systems can be easily reduced to 
the framework of ([T]), Q involving two small parameters g and (3 and characterized by the 
following properties: (i) the linear part is L has large factor ^ before it; (ii) the nonlinearity 
F (U) is independent of g, (3 or depends on g regularly; (iii) the initial data depend on (3 
so that they do not vanish as /? ^ 0; (iv) the solutions are considered on the time interval 
Q < T < where > does not depend on /?. Notice that solutions to ((1]), ([3]) under 
the above conditions exhibit nonlinear effects uniformly with respect to small g, [3 on the 
time interval < r < t*. 

There are important classes of problems which can be readily reduced to the framework 
of H]), ([SI) by a simple rescahng. 

Systems with a small factor before the nonlinearity. Consider a problem of the 
form 

= -iLv + af (v) , v|j^o = h, < a < 1, (110) 

where initial data are bounded uniformly in a. Such problems are reduced to ([1]) by the 
time rescaling t = ta. Note that now g = a and the finite time interval < t < 
corresponds to the long time interval < t < T*/a. 

Systems with small initial data on long time intervals. The equation here is 

dtv = —iLv + fo (v) , v|j^p = aoh, < ao ^ 1, where (HI) 

ioM = it^ (v)+f('"+^' (v) + ..., 

and f^"*^ (v) is a homogeneous polynomial of degree to > 2. After the rescaling v = aoV 
we obtain the following equation with a small nonlinearity 

= h, (112) 

which is of the form of (jllOp with a = a™^^- Note that nonlinearities f in (jllOp which 
are obtained from problems with small initial data and regular nonlinearities fo (v) have a 



dtY = -iLV 



(V) + aof°('"+i) (V) + 
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_ rim) 



special form. Namely, they are almost homogeneous, f (V) = fg 

, m— 1 



term fg™'' (V). Introducing the slow time variable r 
equation of the form ([1]), namely 



(V)- 



with leading 



ton 



we get from the above an 



drV = -- 



rLV 



f (V) + aof (V) + 



it=o 



(113) 



where the nonlinearity does not vanish as ao — > 0. In this case g — a™ ^ and the finite 
time interval < r < t, corresponds to the long time interval < t < with small 

ao <C 1. Note that Corollary [35] for g-dependent nonlinearities can be applied to this case. 
This allows, in particular, to apply results of this paper to Sine-Gordon equation where 



High-frequency carrier waves. Sometimes high spatial frequency of carrier waves 
in the initial wavepackets after a rescaling creates a large parameter ^ at the linear part. 
For example. Nonlinear Schrodinger equation 

d^U = -idlU + ia \U\^ U, C/|^^o = hi iPx) e'^^'=*i^ + /12 (M e'^^'^-^^ + c.c, (114) 

where c.c. stands for complex conjugate of the prior term, and M ^ 1 is a large parameter, 
can be recast in the form ([T]). Indeed, changing variables y = Mx in the above equation we 
obtain 

drU = -i-^dp + ia |C/|' [/, C/|,^o = hi {Piy) e^^''^ + h2 {(iiv) e"^'^^ + c.c, 

where Pi = <^ I, g = ^ <^ I. Note that though the nonlinearity IC/]^ (7 in pi4|) is 
not complex homogeneous, it can be considered as a restriction of a system with a complex 
homogeneous nonlinearity as ([57)1 is a restriction of ([5^ . 

First order hyperbolic equations and systems. Consider now the system (PS)) . 
dlS]) for which the symmetry ([7]) does not hold. The system can be put into the standard 
framework by formally adding two more equations 

Cl C2 
drWl = —dxWi+ Fi{wi,W2) , drW2 ^ —dxW2+ F2{wi,W2) , (115) 

Q Q 

which have only trivial solution wi = W2 = Q not affecting the solutions to the original system 
(|l5l) . (Uni). The extended system has the linear part with two-band dispersion relations 
wi.f (fc) — ciC, |fc|, UJ2X (fc) — C2C \k\ ,C = ±, satisfying evidently ([7]). 



5 Integrated evolution equation 

Using the variation of constants formula we recast the modal evolution equation ([3]) into the 
following equivalent integral form 

U (k, t) = e^^^^^'^^^F (u) (k, r) dr' + e^^C^^h (k) , r > 0. (116) 

Then we factor U (k, r) into the slow variable u (k, t) and the fast oscillatory term as in 
(fn]) . namely 

U (k, r) = e-^f ^(''^u (k, r) , U„,c (k, r) = u„,c (k, r) e'^^^-^''), (117) 
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where (k, t) are the modal cocfScients of u(k, t) as in ([5T|) . Notice that u„_^ (k, r) 
in (jll7[) may depend on g and (jll7[) is just a change of variables and not an assumption. 
Consequently we obtain the following integrated evolution equation for u = u(k, r), r > 0, 

u(k,r)=.F(u)(k,r) + h(k), ^(u) = ^ ^(™)(u'"(k,r)), (118) 

^(™) (u") (k, r) = e^^^'^^F™ (^(e"F^L( )u) (k, r') dr', (119) 

where are defined by and in terms of the susceptibilities x^"^\ ^^^d jr(™) are 
bounded as in the following lemma. 

Lemma 25 (boundness of multilinear operators) defined by 1186]) . U19\} is bounded 

operator from E = C ([0, r*] , L^) into ([0, t*] , L^) satisfying 



^('")(U1...U„,) <T, 

(ui...u„) 



^""ln,., lift 



< 



Proof. Notice that since L (k) is Hermitian, exp |— iL (k) 
inequality 

llu* v||ii < ||u||ii ||v||ii . 

together with ((86]), (fTT9)) we obtain 



(120) 
(121) 

= 1. Using the Young 
(122) 



< sup 



(k,fc) 



(m) 



(ui...u„,)(-,t) 

III |ui(k')|... u™(k(™' (k,fc)) dk'...dk(™-i)dridk< 

(") r||ui(Ti)||^,...||u„(ri)IL. dri<r, 
Jo 

proving ()120|) . Similarly we prove (|12ip by 

\drT^"''^ (ui...u™)(.,r) 
/ / |ui(k')|... u™(k(™) (k,fc)) 



dk' . . . dk(™-i)dk < 



IB • 



The equation (IllSp can be recast as the following abstract equation in a Banach space 

u = J^(u) +h, u,h e F, (123) 
and it readily follows from Lemma [^D that ^(u) has the following properties. 
Lemma 26 The operator T {u) defined by hll^l - UT^) satisfies the Lipschitz condition 

ll^(ui) ~T{n2)\\E < t.Cf ||ui - U2IU (124) 
where Cp < C^rup (4i?)™^^"^ if ||ui||^ , ||u2||^ < 2R, with as ir. 



We also will use the following form of the contraction principle. 
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Lemma 27 (Contraction principle) Consider equation 

x = J'(x) + h, x,hEB, (125) 

where B is a Banach space, T is an operator in B. Suppose that for some constants Rq > 
and < q < 1 we have 

||h|| < i?o, ||^(x)|| <i?o ^/ ||x|| <2i?o, (126) 
||^(xi)-^(x2)|| < <7||xi-X2|| ||xi||,||x2|| <2i?o. (127) 

Then there exists a unique solution x to the equation U25]) such that ||x|| < 2Rq. Let 
1 1 hi II , ||h2|| < Rq then the two corresponding solutions xi,X2 satisfy 

||X1||,||X2|| < 2i?o, ||X1-X2|| <(l-g)"M|hi-h2||. (128) 

Let xi , X2 be the two solutions of correspondingly two equations of the form il25\) with T\ , 
hi and T2, h2. Assume that that T\ (u) satisfies \12b]] . {121^ with a Lipschitz constant 
q <1 and that {{Ti (x) - T2 (x)|| < 6 for ||x|| < 2Rq. Then 

IIX1-X2II < (l-g)-^(5+||hi-h2||). (129) 

Lemma [5^] and the contraction principle as in Lemma [57] imply the following existence 
and uniqueness theorem. 

Theorem 28 Let ||h||^ < R, letr^ < l/Cp where Cp is a constant from Lemma [SS[ Then 
equation U18\) has a solution u E E = C (\0,t^] which satisfies ||u||^ < 2R, and such 
a solution is unique. 

The following existence and uniqueness theorem follows from Theorem 1281 

h 



< R. Then there exists a unique 



Theorem 29 Let ^ satisfy ^ and h e L^ (R'^) , 

solution to the modal evolution equation (0j in the functional space ([0,t*] , i^) . The 
number t* depends on R and . 



Using the inequality (|2ip and applying the inverse Fourier transform we readily obtain 
the existence of an i^— solution of ^ in ([0,t*] ,L°° (M"^)) from the existence of the 
solution of equation ([3]) in C"'^ ( [0, t*] , i^) . The existence of F-solutions in spaces of spatially 
smooth functions can be derived by replacing Lemma [25l with an estimate similar to the one 
in Lemma [50l 

Let us recast now the system (|118l) - (|119p into modal components using the projections 
n„_^ (k) as in pT|). The first step to introduce elementary modal susceptibilities X^™-'*? 
having one-dimensional range in C^'^ and vanishing if one of its arguments Uj belongs to a 
(2 J — l)-dimensional linear subspace in C^"' (j-th null-space of X^™''^ )• For example, in the 
linear case to = 1 when x^^-* acts in C^"^ and is presented in the standard orthonormal basis 
{^nx} i'^ by a 2J X 2 J matrix with elements a^^^, = a^^^^ ^, where index = n,(^ takes 
2 J values, the action of elementary susceptibility x^^^ „/ on a vector v G C^"^ is given by 
the formula x^^^^ n' c'^ ~ '^n \ n' c' ' ®"''C') ^"-C where {e„.^} is the standard orthonormal 

basis in C^-^. Obviously Xnx,7i'X'"^ = n„^f x^^^n„/,(;/v and x''^^^ =J2n.x,n',c' ^n^Cn'X'"^- 
The general definition follows. 

Definition 30 (elementary susceptibilities) Let 

^ = (n, C) G {1, . . . , J}™ X {~1, ir - S", (n, C) e S (130) 
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id X^™^ ^k, [ui (k') , . . . , Um (k*^™))] be m-linear symmetric tensor (susceptibility) as in 



186]) . We introduce elementary susceptibilities x'™^-(k, A: 
tensors defined for almost all k, k by the following formula 



^2.7^'" 



) — > C^'^j as m-linear 



-(k,fc) [ui(k') 



ui(k'),...,u™(k("))] = (131) 



n„,c (k)x('") (k,fc) [(n„,.^, (k')ui (k') , . . . ,n„^ (k(™) (k,^)) u„, (k(™) (k,fc))) 

Then using (I82p and the elementary susceptibilities p31|) we get 



(k, fc) [u, (k') , . . . , u„ (k(™))] . 5:^^^^ (k, k) [u, (ko 



(132) 

Consequently the modal components JF^™''- of the operators jr(™) in (|119p are m-linear 

oscillatory integral operators defined in terms of the elementary susceptibilities (jl32p as 
follows. 



Definition 31 (interaction phase) Using notations from I186\) we introduce for£ — yT-Xj G 
S™ operator 



^(;"^)^.(ui...u„0(k,r) 



Jn. 



exp < I 



|(k,fc) 



(133) 



(k, fc) [ui (k', ri) , . . . , u™ (k(") (k, fc) , n)] d("-i)''MTi, 
wit/i the interaction phase function defined by 

^^k,fcU<^„^,^,,,^(k,fc) 



(134) 



= C^„ (Ck) - C'c^n, (C'k') - ... - C^™^^n„ ((("'k^")) , k(™) = k(") (k, fc) 

Using J^^"'- in (fT33)) we recast J^(™) (u™) in the system (fTT8)) - (ITT9l) as 
^(™) [ui . . . , u™] (k, r) = ^ [ai . . . u„] (k, r) , 

yielding the following system for the modal components ii„_^ (k, t) as in (llip 
"".C (k, r) = Y.^^^^ E,~,H.^ ^) + ^^".^ ' ^) ^ ^ 

6 Wavepacket interaction system 



(135) 



(136) 



The wavepacket preservation property of the nonlinear evolutionary system in any of its 
forms dH), ([3]), (jllSp . (|123p . (|136p is not easy to see directly. It turns out though that 
dynamics of wavepackets is well described by a system in a larger space E"^^ based on the 
original equation (|118p in the space E. We call it wavepacket interaction system, which is 
useful in three ways: (i) the wavepacket preservation is quite easy to see and verify; (ii) it can 
be used to prove the wavepacket preservation for the original nonlinear problem; (iii) it can 
be used to study more subtle properties of the original problem, such as NLS approximation. 
We start with the system (IllSp where h (k) is a multiwavepacket with a given nfc-spectrum 
S — {(k,i,n;), ^ = l,...,iV} as in ([3T|) and fc-spectrum Ks = {k,i, i = \, ...,\Ks\} as in 

(1321). 
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When constructing the wavepacket interaction system it is convenient to have relevant 
functions to be expHcitly locahzed about the A:-spectrum Ks of the initial data. We imple- 
ment that by making up the following cutoff functions based on (j25p 



(k) = * (k, i?k„) = ^ (/3-(^-^) (k - i?k„)) , k„ 6 ifs, « = 1, • • • , li^sl , = ± (137) 
with e as in Definition [T] and /3 > small enough to satisfy 

/3^^^ < ttq, where ttq — ttq (S*) < - min dist {k*i, (Tf,c} ■ (138) 

2 k.iGXs 



In what follows we use notations from (|92|) and 

r=(Zi,...,;„oe {!,... ,ivr, ^= (^',...,^(™)) e{-i,ir, a = (;:^) ga-, (139) 

n=(7ii,...,n™)G{l,...,jr, Ce{-l,ir, (140) 
^ = [n, e , k = (k', . . . , k(")) G M", where S" as in (fnOl) . 

Based on the above we introduce now the wavepacket interaction system 

w,,^ (.) = (., ,?k„,) n„,^ (•) T (E(,,,^,)eA ^''.''') + * (•, ^k, J n„„^ (•) h, ^) e A, 

(141) 

w = (wi,+, wi^_, ...,^N,+,^N-) e -B^^, w;,i; e (^,1?) e A, 

with [■,'dKi) ,^n,-d being as in p37)) . (fTTj) . defined by pl8)) . and the norm in E'^^ 
defined based on ITTl) by the formula 

The index ?9) which takes 2N values labels equations and variables, the right-hand side 
of (|14ip is well-defined for all w G E"^^ and the equality p4ip is understood as equality 
of elements of E"^^ . We also use the following concise form of the wave interaction system 
(fUDl 

^ = T^ (w) + h^, where (142) 

K = («'ii,+n„i,+h,^'i,^„n„,^_h, ...,^'i„,+n„„,+h,i',„,_n„„,_h) ^e^'^. 

The following lemma is analogous to Lemmas [25l [26l 

Lemma 32 Polynomial operator J-<^ (w) is hounded in E'^^ , J-^ (0) = 0, and it satisfies 
Lipschitz condition 

(wi) - (w2)||£;2]v < Ct* ||wi - ^2\\e2« : (143) 

where C depends only on as in I188\}. on the degree of J- and on ||wi||g2jv -f- ||'W'2||^2jv, 
and it does not depend on (3 and g. 

Proof. We consider every operator JF*-™-*^ (w) defined by (|133p and prove its boundedness 



exp 



< 1 and 



and the Lipschitz property as in Lemma [25] using the inequality 
estimates (P5|). Note that the integration in ti yields the factor and consequent 

summation with respect to n,(,i yields (|143p . ■ 

Lemma [32] and the contraction principle as in Lemma [27] yield the following statement. 
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Theorem 33 Let h„ < R. Then there exists Ri > and > such that equation 
Jj^ j[ ) has a solution w G E"^^ which satisfies ||w||^2iv < Ri and such a solution is unique. 

Lemma 34 Every function W;^^ (k, r) corresponding to the solution of ^14^ from E"^^ 
is a wavepacket with nk-pair (k*/, ui) with the degree of regularity which can be any s > 0. 



Proof. Note that according to ()137|) and ()142|) the function 

wz,^(k,r) = ^'(k,79k,,jn„,,^^(k,T), |i^(T)||^i < C, 0<r<T, 

involves the factor 'I';^^ (k) = (^(3^^^^^'' (k — i9k,;)^ where e is as in Dcfinition[TJ Hence, 

n„,^'Wi^^ (k, r) = if n ^ or i?V ^, (144) 
Wi,^(k,r)-*(k,z9k„Jwi,^(k,T),Wi,^(k,r) = 0if |k-i?k,j| > /3l-^ (145) 

and, consequently. Definition [1] for w;^^ is satisfied with Dh — for any s > and C — 

in dsnn- ■ 

Now we would like to show that if h is a multiwavepacket, then the function 



(146) 



is an approximate solution of equation (|123|) (see notation ((92| V To do that we introduce 



^ k^ (147) 



Expanding m-linear operator jr(™) ^ w;^^^ ^ and using notations ([92|l. ([93]) we get 

^^"^ ((E,, -'^^)") = ExeA^ ^'"^ (^a) , where (148) 



wr = wai...wa„, a = (Ai, A„) e A" 



(149) 



The next statement shows that (|146p defines an approximate solution to integrated evolution 
equation (|118p . 

Theorem 35 Let h be a multi-wavepacket with resonance invariant nk-spectrum S with 
regularity degree s, 'w be a solution of ^TJ^ and w (k, t) be defined by ^TJ^. Let 



D (w) = w — (w) — h. 
Then there exists Pq > such that we have the estimate 



D(w) 



<Cq + C(3', ifO<Q<l, P<Po 



Proof. Let 



^ (w) = (i-^^^*,,,^n„,,^j^(w), h =h-^^^*,,,^n„,,^h. 



(150) 



(151) 



(152) 



Summation of p4ip with respect to /, yields 

w = ^ *»,,isn„,,i,JP'(w) + E, ^ *^,,^*^„,,^h. 
Hence, from (I14ip and (|150p we obtain 

D (w) =h- -T- (w) . 



(153) 
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Using (pS)) and ([50]) we consequently obtain 



< CP' if ^ 



< CP' if k„, ^ 



T- (w) = (i - ^^^^ *,,^,n„„,) ^ ^(") (w™) . 

m 



Now, to show (I15ip it is sufficient to prove that 

11^" {^)\\e<C2Q 

Obviously, 
Note that 

Using ([5^ and ( If47p we consequently obtain 

E,.,E(„,,.,,,*(-'^kOn„,, + *^.f, 

with S defined in ((90|) . Let us expand now jr(™) (w™) using (I148|) . According to 
P59p to prove pSSp it is sufficient to prove that for every string A G A™ the 
inequalities hold 

|*oon„,^^(™) (wx)|| < C^3^>for e A, and 

*(-,i9kOn„,^^(™) (w^)|| < C3^?, if (n,k,) e 

We will use (|144p and p45p to obtain the above estimates. According to p35p 
[w,] (k, r) = 5:^^^ [w., . . . w. J (k, r) . 

Note that according to pii)) if = 

wx- = n„_^WA . , if ri = III and i?' = i?. 

Let us introduce notation 

n{tj = {ni,,...,nij, ^ (a) = (n (tj , for A = (/^ i?) £ A™. 

Since 

nri',^n„ ,5' =0, if n 7^ n' or d' ^ d 



(154) 

(155) 
(156) 
(157) 

(158) 
(159) 

fT56l) and 

following 

(160) 
(161) 

(162) 

(163) 

(164) 
(165) 



then (|163p implies 

•^l"c^C ^^^^ ' ' ■ '^^"^ ^ if ^ = [n, C) ^ C (a^ , and, hence. 



^(")[W,](k,r) = 5:^^^^(";;.(,^[w.,...w.J(k,r) 
where we use notation p64p . Note also that 



n„',^^^™l = if n'^"ori?^C, 



(166) 



(167) 
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and, hence, we have nonzero n„'^^^^'"''- (w^) only if 



By Km 



(m) 
,C,?(A) 



(wx) (k,r) = 



cxp < 1 



k,k 



(168) 



(169) 



^^^^-^-^ (k,^) [wA, (k',Ti),...,WA,„ (k(") (k,fc) ,ri)] ^("-D'^fcdri, 
Now we use (|145p and notice that according to the convolution identity in 

|WA, (k', Ti)| • ... • |wA„ (k('»' (k, fc) , n) I = if |k - 

Hence the integral (|169p is nonzero only if (k, fc) belongs to the set 



(170) 



Br, 



< (3^ i = 1, ...,m. 



k-^^t9.k,,,| <m/3i-^}. (171) 



We will prove now that if (n, k*i) ^ S then for small /3 one of the following alternatives 
holds: 



either * (•, Un',^T["l\ (w^j) = 



or (fTO5)) holds and 4>^^^(k,kj 



> c> for (k,k^ e B;3. 

Note then since (t>n ^) smooth then using notation ([M]) we get 

(k, fc) - (k,„ I < CP'~' for (k, fc) e B^, 

(i?i,...,i?„,), k,, ^CJI,'^^!^*'. =C>^m ('?,^) , 
Hence the alternative (|173p holds if 



(172) 
(173) 

(174) 



(175) 



and, consequently, it suffices to prove that either (|172p or (|175p holds. Combining p7ip 
with * (k, z9k„) = for |k - -i?k*i| > (3^^" we find that ^'i^ijj^^™) [w^j] can be non-zero for 

small /? only in a small neighborhood of a point C>!r„i S I'^'l/f out' ^^'^ ^^^^ possible 

only if 

k*, - C^m t) = i?k„, k„ G Ks. (176) 
Let us show that the equality 

Kai^**^''*) =^ (177) 



is impossible for k^^, as in (|176p and rt' ~ n as in (jl67p . keeping in mind that {n, k*j) ^ 5*. It 
follows from ([M)) and (|134p that the equation (|177p has the form of the resonance equation 
(|100p . Since nfc-spectrum S is resonance invariant, in view of Definition 1181 the resonance 
equation p77p may have a solution only if k** = k*i, i = ii, n = ni, with (ni,k,ij e S. 
Since (n,k»j) ^ S that implies (|177p does not have a solution and, hence, (|175p holds when 
(n,k,i) ^ S*. Notice that Theorem 1551 and (|12ip yield bounds 



|wa,|Ib<^i, l|a.WAj|^<C. 
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These bounds combined with Lemma proven below, imply that if (|175p holds then (|16ip 
holds. Now let us turn to (fTBOl) . According to ([T47|) and (fTfO)) the term ^'oon„/,^:?^('") (w^;) 

can be non-zero only if (^>Cm (j^ — k** ^ Ks- Since nfc-spectrum 5* is resonance invariant 

we conclude as above that inequality (|175p holds in this case as well. The fact that the set 

of all (j^ is finite, combined with inequality (jl75p . imply (|173|) for sufficiently small 

/3. Using Lemma [551 as above we derive (|160p . Hence, all terms in the expansion p56p are 
either zero or satisfy (|160p or (|16ip implying consequently (|155p and (|15ip . ■ 
Here is the lemma used in the above proof. 



Lemma 36 Assume that 

vl/,,^,H„,^cxl"|^- (k, fc) [wA, (k', n) , . . . , (k^") (k, fc) , n) 



= 0/or (k,fc) eS^, 



? k,fc 



> tj* > /or (k, fcj ^ B/3, with Bp as in {171^ . 



Then 



*(,/k„)n„,c^;;;).(w3;) 



< 



(178) 



(179) 



4e 



(m) 



n. 



(m) 



■-J " " "-^ CJ* 
Proof. Notice that the oscillatory factor in (|133p equals to 



j 1 1 _E ■ 



exp < L 



Denoting (/)^^^^^ 
Ti we obtain 



(m) 



i(/) (k, fc) 



exp |i0 ^1 



xlj™'' and integrating (|133p by parts with respect to 



*(k,^'k„)n„,,cJ^™; (w,-) (k,r) 



pe 



,0(k,fc)f 



* (k, i?'k,,^ 



(180) 



*(k,i?'k„)- 

S ifl 



(k,fc) 



vi.(k,^'k„)^^^^xr' 



J B 



i0(k,fc) ' 



4"^ (k, k) WA, (k', r) . . . WA„ (k(™) (k, fc) , r) d("-i)'^^ 
xj"^ (k, fc) WA, (k', 0) . . . WA„ (k(™) (k, k) , O) d(™-i)'^fc 

k,fc) dr, fwA, (k')...WA„, fk(") fk,fc))l d(™-i)'^fcdTi, 



where B is the set of k'*^ for which (|17ip holds. The relations (|88|) and (p5| imply 



1""^ tj < Using then P75| .the Leibnitz formula and we obtain ^7^^)- 



The main result of this subsection is the next theorem which, when combined with 
Lemma I34| implies the wavepacket preservation, namely that the solution u„^^ (k, r) of 
P36p is a multi- wavepacket for all r e [0, r»]. 

Theorem 37 Assume that conditions of Theorem \35\ are fulfilled. Let \in,§ (k, t) for n ^ ni 
and Wi^^(k, t) be the solutions to respective systems US 6]) ( [_?.^_?[ ), w be defined by |_?46p . 
Then there exists /3q > such that 



K,^-H„„^w||^ <Ce + C'/3^ /orO</3</3o. 



(181) 
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Proof. Note that u„ ,9 — n„^i5U where u is a solution of (jllSp and, according to Theorem 
EH llull^ < 2i?. Comparing the equations (jllSp and pSOp , which are u = ^(u) + h and 
w = T iyf) + h + D (w), we find that Lemma [23 can be applied. Then we notice that by 
Lemma [26l J- has the Lipschitz constant CpT^, for such u. Taking Cft<, < 1 as in Theorem 
US] we obtain from (fT^ . ■ 

Notice that Theorem[5] is a direct corollary of Theorem [571 and Lemma [M[ The following 
corollary shows that inequality (|181|) and, therefore, Theorems [5] and [3] on preservation of 
wavepackets hold in the case when the coefficients of operator F I U ) in ([3|) , ([86]) regularly 



depends on small g, F ^U^ = F ^U, . 



Corollary 38 (parameter dependent nonlinearity) Assume that conditions of Theo- 
rem \35\ are fulfilled. Consider a perturbed equation mi8\) u (k, t) ^ T (u) (k, r)+jFi (u, g) (k, t)- 
h (k) where operator J-'i (u, g) satisfies the inequality \\J-i{u,g)\\^ < Cg'^ for ||u||^ < 2_R 
with some q, < q < 1. Letv/i -g (k, r) be the solution of ^41\ l- Then l|n„ ,ju — w; < 
Cg'^ + Cp". 

Proof. The statement follows from (|18ip and Lemma [27] ■ 

The following theorem shows that any multi-wavepacket solution to (|118p yields a solu- 
tion to the wavepacket interaction system (|14ip . 

Theorem 39 Let u (k, r) be a solution of U18\) and assume that u (k, r) and h (k) are 
multiwavepackets with nk-spectrum S — {(n;,k^;) ,1 — 1, ...,iV} and the regularity degree s. 
Let also 'I'i,.^ — 'I'i,.iS be defined by ^137^ . Then functions wj ^ (k, t) = ^'j^ ,jn„, ,511 (k, r) 
are a solution to the system |J4J[ ) with h (k) replaced by h' (k, t) satisfying 

h (k) - h' (k, t) < CP", < r < r*. (182) 

Proof. Muhiplying (|118p by we get 

w; ^ = (•, i9k,,, ) n„, (u) (k, r) + * (•, t9k„, ) n„, .^h (k) , w; ^ = ) n„,.^u. 

(183) 

Since u (k, r) is a multiwavepacket with regularity s we have 

||u(.,r)-w'(.,r)||^, <C,r where w' (-, r) = ^^^^ * (•, ,9k„,) u (•, r) . (184) 

Let us recast (I183P in the form 

w;_^ =^'(-,z?k„,)n„,^^.F(w')(k,r) + *(-,i?k„,)n„,,^ [h (k) + h" (k, r)] , (185) 
h" {KT) = [T{n)-:F{w')] (k,r). 

Denoting h (k) + h" (k, r) = h' (k, t) we observe that p55)) has the form of piT|) with h (k) 
replaced by h' (k, r). Inequality ([T^ follows then from PM)1 and pM)) . ■ 



7 Reduction of wavepacket interaction system to a min- 
imal interaction system 

Our goal in this section is to substitute the wavepacket interaction system (I14ip with a 
simpler (minimal) interaction system which describes the evolution of wavepackets with the 
same accuracy. We fix the nfc-spectrum 5* = {(n;,k,i), I = 1,...,N} of the initial mul- 
tiwavepacket and assume everywhere below that it is resonance invariant. The minimal 
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interaction system is built based on operators L and F ^Uj and on S. We want the mini- 
mal interaction system to satisfy the following requirements. Firstly, the approximation of 
solutions of (|14ip by solutions of the minimal interaction system of the order (/i, v) has 
to be of the order g in suitable region of parameters {g, 0) (which is larger for larger /i, v). 
Secondly, the minimal interaction system of the order (/i, v) should be defined by S and by 

the values of L (k) and its derivatives of the order up to /i and by the values x^™^ (k, fc^ 
and its derivatives of order up to v at k*; G Sk- 

The construction of the minimal interaction system consists of the following consecu- 
tive steps: (i) introduction of a time averaged wavepacket interaction system obtained by 
discarding non-resonant terms in the nonlinearity; (ii) reduction of the system for vector 
components v;^^ to an equivalent one for scalar amplitudes vi_^] (iii) change of variables 
k = 'iJk*/ -|- j3ri in the equation for vi^^ resulting in a regular dependence of coefficients on 
small fir]] (iv) substitution of the general dependence on f3ri in the linear part with a cer- 
tain polynomial one of the order /i, and the general dependence on (3r) of coefficients of the 
nonlinearity with a certain trigonometric polynomial of the order (v) substitution of the 
cutoff functions 5* (•j'iJk^iJ from (|14ip . which were preserved up to this step, with 1. 

As a result we obtain a minimal interaction system with weakly universal nonlinearity, 
which in the simplest case, where S is just a single element (k.^, n), is equivalent to the clas- 
sical NLS equation, and in the case when S consists of only two elements (k,t,n), (— k*,n), 
is equivalent to the classical coupled modes system. 



7.1 Time averaged wavepacket interaction system 

Here we modify the wavepacket interaction system ()141|) . substituting its nonlinearity with 
a certain universal or conditionally universal one obtained by the time averaging, and prove 
that this substitution produces a small error of order g. As the first step we recast (|14ip in 
a slightly different form by using expansions (|148l) . ()162p together with p66p and (I167P and 
writing the nonlinearity in the equation ()14ip in the form 

with J^^'^\ ^ as in (|133p and nil \ as in (|164p . Consequently, the wavepacket interaction 
system (I14ip can be written in an equivalent form 

(186) 

The construction of the above mentioned time averaged equation reduces to discarding 
certain terms in the original system p86p . First we introduce the following sets of indices 
related to the resonance equation (|100p and fi™ defined by ([M)) : 

K:^ = {a = (l, C) e A™ : n,^ (d, m, a) = O} , (187) 
and then the time-averaged nonlinearity by 

Note that the nonlinearity ^ (w) can be obtained from by the averaging formula 

()70p where At is defined by formula (|69p with frequencies (f)^ = a;„^. (k*i^. ). Consequently, 
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the desired equation with time-averaged nonhnearity is 

V,,^ = * (•, l9k,,J ^av.n,,^ (v) + (•, ^9k,,J n„,,^h, I = 1, ...N, ^ = ±, (189) 

which similarly to (|142p we recast concisely as 

v = ^av.*(v) + h^. (190) 
The following lemma is analogous to Lemmas [321 1211 

Lemma 40 Operator J^av.* (v) is bounded for bounded v G E"^^ , ^av.* (0) = 0. Polyno- 
mial operator J^av.* (v) satisfies the Lipschitz condition 



* (vi) * {^2)\\e2N < Ct^ II Vl V2||^2]\ 



(191) 



where C depends only on a in \89(l . on the power of T and on ||vi||^2jv + ||v2||£2jv, and, 
in particular, it does not depend on (3. 

From Lemma HOI and the contraction principle we obtain the following Theorem similarly 
to Theorem [551 



Theorem 41 Let 



E2N 

2N 



< R. Then there exists i?i > and > such that equation 



il90]) has a solution v G -E^^ satisfying \\'v\\^2n < -Ri, and such a solution is unique. 

Theorem 42 Let vi^^ (k, t) be solution of U89\) and w;^,? (k, r) be the solution of ^Jl^. 
Then the v;^^ (k, r) is a wavepacket satisfying ^TJ^, (lj5^ with w replaced byv. In addition 
to that, there exists (3q > such that 



||v/,^ -W/,^||£<Ce, l^l,...,N; ^ = ±, forO<g<l, 0<P<Po- 



(192) 



Proof. Formula p44p . (I145P for v/^^ (k, r) follow from (|189p . We note that w is an 
approximate solution of (|189p . namely we have an estimate for Dav (w) = w — .?"av,* — hi^ 
which is similar to 



Dav (w; 



w — 



av,* 



<Cg, if < (? < l,/3 < 



(193) 



The proof of (|193p is similar to the proof of (|155p with minor simplifications thanks to the 
absence of terms with 'i'oo- Using (|193p we apply Lemma \Tl\ and obtain p92p . ■ 



7.2 Averaged system for scalar amplitudes 

Now we recast (|189p in the form of an equivalent system of scalar equations for amplitudes 
vi^^ = v\ of solutions va, defined based on (jlip . namely 

VA, (k) = VI/ (k, C^'^k,,,) n„,,^(o (k) VA, (k) = (k) g„,_^(o (k) . (194) 

Note that according to p45p support of ^(o is localized near Ck*i, ,and we can assume that 
Srii c''' depend smoothly on k near this point. Multiplying (|189p by grij,^; (k) (with the 
standard scalar product in C^^) and using p94p we obtain the following system of scalar 
amplitude equations 

= ^'(•,i?k„J/av,„,,^ (w) + «'(•, 'i?k*ij/i„,,i,, l = \,...,N, ?9 = ±, where (195) 
K.:0 - M ■ n„„^h, /av,„,,^ {v) - E„,ean. ^A-eAj'^,, Gm ^"^^^ ' ^^^^^ 
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According to (|169p the m-linear operators in the above equation are given by 



/"^-(^j) (k,r) 



Q^^l^ (k, fc) = g„,, (k) . (k, fc) [g., (k') , . . . , (k(™) (k, fc) ) 
The concise form for the system (jl95p of scalar equations for ampUtudes is 



(197) 
(198) 

(199) 



where the components vi^f) of v belong to the space i?sc of scalar functions with the norm 
defined by (fT7)) . (|18p applied to scalar functions. Note that Q^™^ - ( k, fc ) can be extended in 

an arbitrary way as bounded functions for arguments k, k where (I17ip is not satisfied, for 
example the extension can be zero, the extension does not affect solutions of (|195[) because 
this equation involves factors ^' (^i^k^ij and (|145p holds. 

Lemma 43 Operator is bounded for bounded v G E'^^ and (0) — 0. Polynomial 
operator /if (u) satisfies the Lipschitz condition 

II/* (Vl) - f^S, {v2)\\e2N < Ct^ \\VI - W2||£2« 

where C depends only on as in i88\) . on the order of as a polynomial and on \\vi \\e2n + 
\W2\\e2n, and it does not depend on (3. 

From Lemma HOI and the contraction principle we obtain the following Theorem similarly 
to Theorem [33l 



Theorem 44 Let 



< R. Then there exists i?i > and > such that the \199\] 



has a solution v G E^^ satisfying \\v\\e2n < Ri, and such a solution is unique. 

7.3 Rescaled amplitude equations 

According to (|145p amplitudes vi^^ (C^*/ + rj) are localized about the point rj — 0, and to 
study its behavior in a vicinity of rj = we introduce a group of dilation operators 

(Bpv) (rj) = (3% iPv) , > 0, (200) 
which preserve the L^-norm and commute with the convolution, i.e. 

\\Bpv\\Li - II«IIli , Bpv * Bpw = Bp{v*w). (201) 
We introduce then a rescaled and shifted version of initial data /i„, in (|196p by the formula 

i?„,,^ (k) = Bph„^^^) (k + i^Ki) , K.^-ff (k) = r'^Hn^s (k - ^Ki)) , (202) 
where Bp is defined by (|200p . |k — i?k*/| < /3^~^, and new variables 

r7, = p-' (k - ^Ki) , Z = 1, iV, 77 = (J7i, . . . , r,^) . (203) 

In this and the following sections we assume that ifn,,^ {/3, rj) are defined for all rj E R"^, in- 
cluding \ri\ > f3~'^. Though (|195p involves with a cutoff factor, namely 5* (k, i?k,ij /i„,_^ (k) = 



assume that 



(k) as in we will later use iPiV) defined for all t;, and 



(l-vl/(/3^T,))ff„,,^ if3,rj) 



< CP' 



(204) 
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where (i) ^' {13'^ rj) = (ry, 0, /3~'^) is as in ([^5]) . (pS)) : (ii) e and s are the same as in Definition 
m (iii) condition ((204)) is consistent with (HH) and pO)) . 

For a solution {); ,5 (k, r) of (|195p using p45p we introduce the foUowing functions 

zi,^ (77, r) = /J'^Wi,^ (i9k,i + /3j7, r) , z,,^ [r], r) = * (/3^J7) i,,^ (r,, r) , r; e R'', (205) 

which satisfy a rescaled version of (|195p provided below. Note that since ^n, = X E A™ ^ 
and the nfc-spectrum S is resonance invariant we have ^A^ = (^'•^''k*;. = Ck*; = ^?k*/. 

Since k, /c satisfy the convolution identity (|57|) the variables r/, f/ defined by (|203p satisfy 
similar identity as well, namely 

V = T."^,V^'\ r,M(k,r7)=T,-^""\w. (206) 
^ — ^i—i ^ — ^'i— 1 

Change of variables (|203p in the integral operator /av.ni.i? defined by (|197p yields the following 
amplitude system for zi^^ which is equivalent to p95p : 

i,,^ (r?) = * iP'rj) /av,.,.^,^ (^) (r?) + (/3'J7) ^n,.>3 (^) , ^ = 1, -N, = ±. (207) 
According to P?7)) . (jTM)) and PM|1 

* (k, z?k„, , Z?!-^) = VP (/3^^) , /,,_„^^^^^ (f) ^ ^^^^^ /i";^,,,,^ (z) , (208) 



Qlri,l(A-) i^^*' + ^* + ^'^) liZi (''^'O d^""'^ Vn. 



Note that the condition (llTip on the domain of integration takes in the new variables the 
form 

" ' < i = 1, ...,TO and |?7| < m/3"^ (210) 



Finally, we rewrite amplitude system (j207p in the concise form 



z 



= * (/3^-) + * iP'-) Hp, z e Ei^ . (211) 



Let us show now that (|21ip is of the form of (jllSp with 2 J-component vector u substituted 
with 2A'^-component vector z, the matrix L (k) substituted with a diagonal matrix L with 
entries "diVm (i^k*; + Prj). For that we introduce S-averaged tensor defined on z G 

^2Nm ^Yie formula 

iPv, PV, z) = T.XeA-, (A) (^^*' + ^* + f^"^) (212) 

which depends on S through A™^ and acts from C^^™ into C^^. Note that z\. and Q''"^^ ^ 
are scalar factors, za^s a scalar projection in C^^ onto a line along A^-th eigenvector of 
L. Hence, the right-hand side of (|212p is a sum of elementary susceptibilities obtained from 
qI™^ as in p32p and p07p has the form of p36p . Note that non-zero terms in ()212p contain 
products z\. which satisfy (jlOOp . Therefore, if /? = and 5 is resonance invariant qI™"* has 
the form of weakly universal nonlinearity; if S is universally resonance invariant then qI™^ 
has the form of a universal nonlinearity as in (|65p . 
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7.4 Amplitude system with polynomial dispersion relations 

Now we introduce amplitude system with polynomial dispersion which is similar to (|207p 
and provides (i) sufficiently accurate approximation to (|207p : (ii) standard polynomial de- 
pendence of coefficients on t], fj in the sense clarified below. The amplitude system has the 
form 

■i« = » (/J'l) /t;,"' (S) + * (fl'i?) A..„«, ! = 1, -N: i) = ±, (213) 

where 5* (/3*^i7) are cutoff- factors defined in (I208p . (I137P and approximations for 

f are defined below. The indices u = 1,2, = 0, 1 determine the order of ap- 

proximation: (i) fi determines the order of approximation of the dispersion relation by a 
polynomial of the degree /i; (ii) v determines the order of approximation of the susceptibility 
coefficients ()198p by a trigonometric polynomial of the degree ly. As before, we recast ()213p 
in a concise form 

u^'i'fjf'-'''''^ (u) + '^f}H. (215) 

where (tj) = {P'^rf). Finally, we eliminate in (|213p the cutoff factor ^E* {P'^rf) by setting 
^ {P'^'n) = ^' (0) = 1, and introduce the amplitude system with weakly universal nonlinearity 
and polynomial dispersion without cutoff 

iv) - /l^,;^ (u) iv) + K.a 1, -N, I? = ±, (216) 

which can be written in the form of (j215p with ^E*^ = 1. 

Let us turn now to the construction of the approximations. For every nfc-pair (k,;,rti) 
we introduce the Taylor polynomials of order /i of the dispersion relation lj„j (k,; + Prj): 

72 {Ki,ni, /3t]) = 7i {Ki,ni,f3r]) + — {t], uj[[^ (Ki) t]) , 
and similarly 73 for /i = 3. Obviously we have the inequality (see p7ip ) 

\cUr.,{Ki+Pv)~l^{Ki,ni,Pr])\ <C/3(^+i)(i-^i), (k, fc) e Bp. (217) 

The phase function (f'n ^ 1 ^ = (j^i > defined by (|134p . is approximated then by a 

polynomial phase function 

cj,[^^]^.^{CKi,t,f3ri,f3f]) = (218) 
C7^ iKi,m,Pri) - C'7^ {Ki, , n', /3r,') - ... - (('"^7^ (k*i„ , n^"), /3r,('")) . 



Note that since ^ = ^Aj with A e A™ defined by ^TE7]i . the equation (fTUDl) is fulfilled. 
Hence, c/)'^''^ - ^i?k*/. A:*, 0, 0^ = and the function (f)^ depends linearly on 77,7/ and 



^ ^ ^ is quadratic, namely 



6^ 



{^Ki,l,Pv,f3fl) = (^k,;,E,77,,7) , (219) 

,z, £, /3r7, /3r/) = /J^^^ (^?k,z, £, r?, 7?) + /3'0';'^_^- (^?k,z, £, 77, 77) . (220) 
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In the case ^ — 2 the polynomial phase function involves two parameters g^, q^'- 

(221) 

_ Q_ 



£-1 = e'2 = ^^2; < ^1 < < ^2 < 00, (222) 



where and ^2 '^^Y be large or small depending on relation between q and /3. Sometimes 
it is convenient to consider and Q2 as independent parameters. If /i = 1 wc formally set 
Q^ — 00, e7 ~ ^- P7ip holds we have the estimate 



< CtJ- , ^ = 1,2. 



(223) 

To ensure that the approximation error is small for given /i we assume that q and /3 satisfy 

(^+i)(i_,) 

(? ^ 0, /3 ^ 0, y 0. (224) 

Now we approximate the dependence of Q*'™''- (i9k*/ + /3?7, fc, + /^Ty) on r;, 77 given by p98p 
by trigonometric polynomials. Zero order approximation with = is given by 

^1"'? l'^^*' + £ + /Jt?) = (i?k,i , E) . (225) 

To define the first order approximation we modify the standard Taylor expansion using 
trigonometric polynomials instead of algebraic ones. Taking the first derivative with respect 
to /3 at /3 = 



Qj;^^|'-(^?k,i,r7,£, 



d 

d/3 



[3=0 



which obviously is a linear function with respect to 77, 77, we express then r] in terms of fj 
using (|206l) : 

Then the first order approximation is 

m 

gj;"^'l^ (dKi + (3rj, K + /Jt?) = (t?k,,, £) + ^ (7j^"j)| (79k,z, £) • sin/Sr?^^"), 

where sinry^-') = ( sinT^^-''' , sinr;!^'''' ) . An advantage of this approximation is that the 



multiplication by sin is a bounded operator which equals the Fourier transform of a finite- 
difference operator whereas the multiplication by 77^"'^ corresponds to the partial derivative 
and is unbounded. Since the original nonlinearity does not involve unbounded operators 
the use of bounded operators is natural and convenient. In fact, it is well known that 
the presence of the derivatives in the nonlinearity of NLS-type equations causes well known 
technical difficulties, see |14j . In our approach the approximating equation provides the same 
accuracy and its nonlinearity involves only bounded finite-difference operators bypassing 
those difficulties altogether. 
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According to Condition [12] the susceptibility is smooth and if (|210p holds we have the 
following inequality 



< C/3 



(l/+l)(l-€l) 



(226) 



We introduce components ,f '^'^ of the weakly universal nonlinearity /('^•'') by the formula 



^77' + ...+77("^)— 77 

m 



As before, we establish standard properties of the operator /(^''') defined by the above 
formula. 

Lemma 45 Operator 'i'pf'^^''^^ is bounded for bounded u G , /,p (0) = 0. Polynomial 
operator 'i'pf^^'"^ satisfies the Lipschitz condition 



(228) 



where C depends only on a in I188\) . on the power of J- and on \\ui\\^2n + ||'U2||£;2jv ■ In 
particular, it does not depend on (3 > and on < < 00, < g2 < 00. 

From Lemma HOI s^nd the contraction principle we obtain the following Theorem com- 
pletely similar to Theorem [33l 



Theorem 46 Let 



< R. Then there exists i?i > and > such that equation 



-B2« 

\190\) has a solution z £ E^^ satisfying ||z||^2jv < Ri- Such a solution is unique and 
z/,^ (k, r) = 1/ |k| > 

Theorem 47 Let iii^^ (k, t) be a solution to i215]) and zi^^ (k, r) be the solution of i211]) . 
Then the following inequality holds: 

\\ui,i, - zi^4e^^ < C/3(^+i)(i-') + Ce-i/3(^+i)(^'^\ / = 1, iV; ^ ^ ±, (229) 

for all < g < 1 and < /3 < /3q, where e is the same as in Definition[l\ /3q is sufficiently 
small. 

Proof. To obtain (|229p we note that ui^.g is an approximate solution of (|21ip . namely 
u - '^pf'^^'"'^ (u) -h^^D where D is small. 



To estimate 



D 



observe that integrals involving u have the integration domain as in (jl7ip . 



Hence, using (|226p and (|223p we obtain 



D 



B2A 



< C(3 



(M+l)(l-e) _,_ C'p-l/3(A'+l)(l-<^) 



and applying Lemma [27] we get (|229p . 
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7.5 Decay of solutions and elimination of cutoff factors 



In this subsection we show how to remove the cutoff function in (|213p and to obtain the 
averaged interaction system with a weakly universal nonlincarity. If /i = 1, = and the 
nfc-spectrum S is resonance-invariant, the amplitude system coincides with the system (1621) 
with a weakly universal nonlinearity. For fi > 1 or v > the amplitude system involves 
additional terms. In particular, ii ^ = 2, v = and S — {(k*,n)} is just a single element 
then the linear part has the second order and the nonlinearity is universal, and amplitude 
system turns into the classical NLS system: 



drU(; = C-72 {K,n, -iCPVr-n) + bQulu.Q, (0) ^ He, C 



This system is equivalent to ([5T|) when i?_ = H^, b- — b'^, — u'^. When v > the 
nonlinearity involves additional terms with finite difference operators. 

The possibility to remove cutoff functions is based on the fast decay of u (k) as |k| oo 
which is equivalent to high smoothness of u (r). The factor "if (3 can be replaced by 1 with a 
small error when data H (k) decay sufficiently fast. To to describe the decay we introduce 
weighted Banach spaces of scalar functions H (k) described as follows. 

Definition 48 (weight function) For a > we call a positive function ip{r), r > 0, 
a weight function from class W (a) if it satisfies the following conditions: (i) 1/^(0) > 0, 
t/i (ri) > ijj (r2) for ri > r2 > 0; (ii) ifj (ri + 7-2) < i) (ri) + ijj (r2) + C where C does not 
depend on ri,r2 (ip is sublinear); (Hi) ip (r) — alnr > C > for all r > (ip (r) is 
superlogarithmic ) . 

We introduce {ip) as a space of scalar functions H (k) , k e M'^ with the norm 



H 



dk. 



(230) 



For vector-functions we use the same formula with Euclidean norm | ■ | . In the simplest case 
of '0 (r) = aln(l + r) we have ip G W (a) and obtain (ip) = L^'° with the norm ([19]). 
If the weight function belongs to W (a) for all a the space (ip) consists of the Fourier 
transforms of infinitely smooth functions. The following Lemma shows that (ip) is closed 
with respect to the convolution. 



Lemma 49 Let Hi,H2 G {ijj) and 



^3(k) 



Hi (k- k')^2 (k- k') dk'. 



Then H3 (k) < C Hi (k) Hi (k) 



(231) 



Proof. Using Definition 25] (ii) we obtain 



p'A(lkl) 



i?3(k) 



< / e'^d'^l) 



c / gV(|k'|)>(|k-k'|) 



ffi(k-k') ff2(k') 
i?l(k-k') i?2(k') 



dk' < 
dk'. 



Applying Young's inequality (|122p we obtain 



Hsik) 



dk < e*^ / e^d^^l) 



Hi{k) 



dk' / e^d'^l) 



i?2(k) 



dk' 



implying (|23T]l . 
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Let us introduce the norm in the space Esc ("0) by the formula (|17p 



Hi;-) 



Hi;-) 



= sup 

C([0,r.],Li(^)) o<r<r 



.^A(|k|) 



HiKr) 



dk. 



(232) 



Using (|23ip instead of we obtain as in Lemma [55] the following statement. 

Lemma 50 Operator ^pf'^-"'' in i215\) is bounded for bounded u G Ej^ (V^), / (0) = 
and satisfies Lipschitz condition 



^/s/^^'"' (^i) - ^/^/^""^ ("2) 



< Ct^, \\ui - U2I 



(233) 



where C depends only on C-^ a in i88\). on the power of polynomial /(^^'^^ and on 11^1 ||£2jv(^)H 
||£;2jv(^) and does not depend on [3>Q and on Q < < 00, < ^ 00 . 

From Lemma |40] and the contraction principle we obtain the following Theorem com- 
pletely similar to Theorem 1331 



Theorem 51 Let 



H 



p2N 



< R. Then there exists Ri > and r, > such that equation 



11215]} has a solution u € E^^ iip) which satisfies \\u\\i^2N(^^-^ < Ri, and such a solution is 
unique. 

The following lemma shows that ^' can be replaced by one with a small error. 



Lemma 52 Let H 
[204\ l holds. 

Proof. We have 



<C,'4;eWia), as in (2^. If s > 0, e > and f < a then 



(1 - * (/3^,7)) Hirj) 



\v\>p- 



Hin) 



dr] = 



7,|>/3- 



,i'{\v\) 



Hiv) 



dr] 

(234) 



< / e-^(^'') e^^^'^^^Hiv) drj < /3'eHP-')'/^-H('"') 

\r,\>l3-' 

According to Definition [48] (iii) 

In {p") s/e - (r") < aln - < C 

and we obtain (|204|) from (l234l) . ■ 



H 



Theorem 53 Let 



H 



^ < R where weight function ijj belongs to W (a) and let ^ < a. 

Let u and uq be solutions to respectively the minimal equation with cutoff factor and without 
cutoff factor respectively. Then there exists Cg and /3g such that 

N- i/o||£2«(^) < Csl3\ 0<P<Po- (235) 
Proof. We show that u is an approximate solution to uq = f^^^^'^^ iuo) + H. Namely, 
u = *^/(''''') iu) + ^pH = /(^'"^ iu) + H + D, i'fp - 1) /(^''^) (u) + i'i'p - 1) H. 

According to Lemma l49l if u G E^^ (V') then /(^''^) (u) e E^^ itp). Applying Lemma [52] we 
obtain 

D < CP\ < /? < /3„. (236) 

Lemma [27l combined with (|236| yields (|235|) . ■ 

Now we give the theorem on approximation by solutions of a minimal system without 
cutoff. 
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Theorem 54 Let Hi^q (k), / — 1, N he functions hounded in (ip) where ip helongs to 

W (a), let f < a. Let hi^Q (k) be defined by \202\) and ^hi^t^ (k) -^hi^i^ (k)g,i,,c (k)- Let 
u (k, r) be a solution of equation mi8\} with multiwavepacket initial data of the form i3S\) . 
Let ui^{) (k, r) he a solution to the system with a weakly universal nonlinearity i216]) with 
initial data ui^g (k, 0) ~ Hijs (k) and 

Umin (k, r) = ^i^^ P^'^ui^^ {l3^^ (k - Ck*z, ) , r) gn,,^ (k) . 

Then 

\\u - u„,i„||^ < + C/jf-^+Dd-') + Ce-i/3(^+^)(i-'^ + Cg. (237) 

Proof. We take u — J2§J2fLi^i,-6 ^^'^ estimate ||u (k, r) — u„iin (k, r) ||^ applying 
subsequently Theorems [371 HH formulas (|194p and (|205p , Theorem 27] and finally Theorem 
[S3] to obtain inequahty (|237p . ■ 

Note that Theorem [7] is a direct corollary of Theorem [54l 

Remark 55 A^ote i/iat i216\) is the Fourier integral version of the following system of equa- 
tions based on weakly universal nonlinearity and slightly more general than i6^) 

1 z ( \ 

drui^n = — w^, (k,i) • VxUi,§ + - — • w"^ (k,i) VrUi^n + fll^^Y '5^) ' (238) 

Wiislr=o = -^'iJ' where dm (r) = (r + e^) - Uj (r - e,;) 

where g^, Qi are as in S222\) and et is i-th standard ort in M.^ . In the case when 15^)) holds 
1/^2 is hounded or small and the dependence on the coefficient 1/^2 is regular for small g 
and [3 and u^j (k, t) may be looked at as a shape function. When gi = g and l/g2 is 
substituted by zero we obtain an equation exactly of the form \6l3jl . 

When = 0, /i = 1 and the nfc-spectrum S is universally resonance invariant as in 
Definition llSl the nonlinearities /^^'^^g are universal of the form (|65p . When the n/c-spectrum 
S is resonance invariant but not universally resonance invariant, the nonlinearities are weakly 
universal, but may be not universal, that allows, in particular, for the second and the third 
harmonic generation. 
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